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PREFACE 


The one day colloquium "Magnetic Field - Plasma Interaction on the Sun" 
was organised at Kodaikanal on January 23, 1984, where the main solar research 
facilities of the Indian Institute of Astrophysics exist and most important, 
the place where John Evershed F.R.S. lived and worked for seventeen years. It 
was on January 5, 1909, seventy five years ago, that he discovered the radial 
motion of the solar plasma in the penumbral regions of sunspots, which was the 
first visible demonstration of the interaction of plasma with magnetic field. 
Originally, the title of this colloquium was planned to be "Evershed Effect", 
but subsequently it was changed to the present, one to provide a wider scope so 
as to cover a range of allied phenomena on the Sun. 

The entire financial support to organise this colloquium and to bring out 
these proceedings as a special issue of the Kodaikanal Observatory Bulletin was 
provided by Prof, J.C.Bhattacharyya, Director, Indian Institute of Astrophysics 
and we would like to express our gratitude to him. We thank the authors as well 
as the participants who have contributed to the colloquium. Dr. S.P.Bagare with 
the staff at the Kodaikanal mountain contributed in many ways to the local arrang¬ 
ements. He also helped us in editing the manuscripts. Miss K.Padrnavathi did 
the entire, typing work and Mr. R.Krishnamoorthy made the camera ready copies 
of the manuscripts and was also personally responsible for the printing and prepa¬ 
ration of the proceedings in their final form. We acknowledge with pleasure 
.the valuable and devoted services of all these individuals, which have made.this 
mission a success. ' • . • 


Bangalore, December 1984 


K,R. Sivaraman 
P,K. Raju 
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MESSAGE 


I am very happy to learn about the Colloquium being 
arranged at Kodaikanal on the occasion of the 75th 
anniversary of the discovery of the Evershed effect,, The 
contributions of John Evershed to the development of the 
Kodaikanal Observatory are overwhelming. He equipped the 
Solar Physics Laboratory with fine instruments, almost 
entirely built with his own hands. The high resolution 
spectrobeliograph set up by him at Kodaikanal was instru¬ 
mental in arriving at answers to some vital questions 
about the gases forming sunspots* His discovery opened up 
a new area in our understanding of the physics of the sun. 

Over these past seventyflve years we have moved into 
many new areas of this fascinating subject. More details 
of the interactions of solar plasma with the magnetic 
field have come to light. Deeper thinking aimed at unveiling 
the energetic processes happening in the solar atmosphere 
are now called for. I hope that the discussions will open 
new pathways to future investigations as the discovery by 
Evershed has done 75 years ago. 

X wish the Colloquium all success. 

(M.G.K. Menon) 

Chairman, Governing Council 
Indian Institute of Astrophysics 
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WELCOMING ADDRESS 


J.C. Bhattacharyya 
Indian Institute of Astrophysics 
Bangalore 560034 
INDIA 

Friends and Colleagues, 

I welcome you all to Kodaikanal and to this Observatory) I wish to express my gratitude 
to all the individuals who have come inspite of their heavy pressure of work? my special thanks 
to our friends from abroad, who have been generous enough to accommodate this visit in their 
very tight schedule. 

The purpose of today's gathering is to discuss some aspects of study of the physics of 
the sun; the occasion is to commemorate a discovery which was achieved in this Observatory three 
quarters of a century ago. On your way up here, you might have noticed a plaque on the wall 
of a building indicating the event. I shall try to give a little history of the humarl efforts 
which culminated in the unfolding of one of nature's mysteries. 

Prof. M.G.K.Menon, Chairman- of the Governing council of this Institute was very keen 
about arranging this meeting. He had planned to attend the meeting but unfortunately, the extreme 
load of his responsibilities have prevented him f'rom coming here today. He has, however, sent 
his message on this occasion. 

Let me start from the very beginning. In the last decade of the eighteenth century, 
a small observatory was started by the British East India Company in the outskirts of the settle¬ 
ment around Fort St. George which has now grown to be the cosmopolitan city of Madras. They 
had, among other instruments, a transit instrument which was mounted on a granite pillar. On 
the four sides of the pillar, the founders decided to leave a message to the posterity. Inscribed 
in four languages, Latin, Tamil, Telugu and Urdu, known at that time as Hindustani, the pillar 
still proclaims the historic event. Translated into English the message reads: "Posterity may 
be Informed a thousand years hence, where the seeds of Mathematical sciences were planted by 
the British liberality in Asia". 

The dreams of the founders were not in vain, as shown by the course of events which follow¬ 
ed. Under the inspired direction of some devoted astronomers, the Madras Observatory indeed 
blossomed into an active centre of astronomical research. Names of Taylor and Goldlngham may 
be mentioned, of those who undertook precise determination of positions of several thousand 
stars. Seven decades later, the helm of affairs passed on to an astronomer whose name is inde¬ 
libly marked in astronomical literature; Norman Pogson, who took over as the Chief Astronomer 
of Madras Observatory. He introduced a fresh wave of enthusiasm and encouraged local, scientists 
in the investigation of the heavens. Name of C.Raghunathachari could be mentioned as one among 
the first Indian astronomers of modern era. Pogson was keen on the study of asteroids. To 
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the list of asteroids, five names, Asia, Sappho, Sylvia, Vera and Camilla were, added by his 
efforts from the Madras Observatory, He put up a strong case for modernization of observing 
equipment, and to shift the observatory to a more suitable location on the hills. 

To cut a long story short, things began to move only after Pogson'o death. A location 
in the Pa Ini hills was selected and the work on construction began in 1899. In 1901, the Kodal- 
kanal Observatory came into existence. The emphasis was now on Astrophysics, the new astronomy. 
Detailed studies of solar spectra had started giving totally new types of information. Only 
a few decades ago, something which was considered Impossible to know about the celestial objects 
had suddenly become available in unimaginable detail. Spectroscopy had become the most powerful 
tool In Astrophysics. 

Evershed came into this field at this juncture. In 1905, he left his village home in 
Surrey in the United Kingdom, where he had built a small but well equipped observatory. He 

packed up some of the components and brought them with him to Kodaikanal. On his way to India, 

he travel led via the United States and Japan where he picked up some more components and came 
fully prepared to start a modern observatory. He had discussed his plans with stalwarts like 
Pickering, Barnard, Frost and St. John and hence he went straight ahead with his task at Kodai¬ 
kanal. The spectroheliograph in the next building that you can see today has basically remained 
in its original configuration designed by Evershed. 

This was not the first visit by Evershed to India; he had earlier been with the British 
team at Talni near Ratnagiri in January 1898 when he first obtained the near Ultraviolet spectra 
of the Sun. On joining Kodaikanal, the first problem he undertook for investigation was to 
explain the peculiarities of sunspot spectra. In the course of that work, fortune smiled upon 
him. On January 5, 1909, after seven days of heavy rain, the sun rose in an unusually clear 

sky with excellent seeing. Evershed could obtain a series of high resolution spectra across 

a spot group. The original plate is preserved in this Observatory; we have placed a magnified 
image of the plate in the spectroheliograph room. By employing the method of photographic sub¬ 
traction, largely developed by him, he could get clear indications of systematic Doppler shifts 
in the spot regions. One of the nagging questions about the dynamical structure of sunspots 
began to clear up. 

That was the beginning of a series of investigations, which unfolded the system of complex 
flow around sunspot regions. Since then it has been a long journey along the path of discoveries 
in many observatories, connecting the plasma interactions with magnetic field, small scale velo¬ 
city fields, large supergranular motions, and many other features displaying an infinite variety 
of plasma dynamics in the solar envelope. 

The problem has not been completely solved today. No scientific discoveries clear a 
problem completely. They give deeper insights into nature's processes and we come across fresh 
problems of a more refined variety. There are several aspects of the phenomena which remain 
as enigmatic as ever. We have chosen today to have a general reconnaisance of the progress 
so far made in this field and of the problems we face today. We hope to obtain some ideas 
regarding the immediate steps needed and therefore to readjust the course of future investigations 
in this field. We shall now begin with the scientific papers. 


xiv 
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THE ROTATION OF THE SUN FROM MOUNT WILSON SUNSPOT MEASUREMENTS 

■k 

Robert Howard 

Mount WiIson and Las Campanas Observatories 
Carnegie Institution of Washington 

and 

Peter A. GiIman 

High Altitude Observatory ** 

National Center for Atmospheric Research 


1. Introduction 

We have recently completed the measure¬ 
ment and reduction of 62 years of white-light 
solar images taken at the Mount Wilson Observa¬ 
tory. So far the data have been analyzed 

for differential rotation and time variations 
of this quantity, among other parameters. 
The data, the measurement and reduction techni¬ 
ques, and the results obtained so far are 

described in two publications: (Howard, Gilman, 
and Gilman, 1984? Gilman and Howard, 1984). 
This is a brief review of that work. 

The interval covered by these data 
is 1921-1982. The measurements were carried 
out in an interval of two years by Pamela 
Gilman. The homogeneity in the reduction 

and in the observations themselves are two 
advantages of these measured spot positions 

over other data. A further advantage of 

this analysis - and a characteristic that 

distinguishes it from nearly all previous 

WO rk - is that measurements were made of 

individual sunspots, not just sunspot groups. 

# 

Now at National Solar Observatory, Tucson, 
Arizona, U.S.A. Operated by the Association 
of Universities for Research in Astronomy, Ino. 
under contract with the National Science Foundation. 

** 

The National Center far Atmospheric Research 
I s sponsored by the National Science Foundation. 


2. Data 

The white-light solar images have been 
obtained with the same focal-plane shutter 
mechanism at the 60-foot tower telescope 
at Mount Wilson throughout the interval. 
The 18.3 m focal length lens now In use was 
originally installed about 1916. The coelostat 
mirrors were replaced twice during the Interval, 

The image diameter is 170 mm, Pole 
markers on the images designate the axis 
of the solar rotation. The focal-plane shutter 
gives an exposure time of about 1/1000 second 
for the f/60 beam. The observing log contains 
the time and data of each observations, the 
focus setting, temperature and other informa¬ 
tion. 

3. The Measurements 

An electronic measuring pad placed 
over a light box was used for position measures 
on the plates in two dimensions. The accuracy 
of such measurements is about 0.1 mm, which 
is about one arc sec on the plates. 

Each plate was positioned on the tablet 
with the rotation axis of the sun aligned 
with one axis of the tablet. The limb was 
defined by 8 equally spaced measurements. 
Then each spot umbra (penumbrae were not 
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mesured) was measured with two position measur¬ 
es of the cross-hair, This enabled us to 
calculate the position and approximate area 
of each spot. 

Altogether nearly 15,000 usable images 
were obtained in the 62-year interval. Nearly 
335,000 daily sunspots were measured on these 
plates. 

On ly sunspot umbrae were measured 
in this program because we felt that they 
represented the fundamental positions and 
areas of the spots. 

4. The Analysis 


The spot group and individual spot 

rotation rates discussed in this paper are 
sidereal rates, 

5, Results 

The spot group rotation rate deter¬ 
mined for the whole data set is: 

w = 14.393 (±0.010) - 2.946 (±0.090) X 

? -1 

Sin B deg day , 

where B is the heliooraphic latitude. For 

all spots we obtained 

to = 14.522 (±0.004) - 2.840 (±0.043) X 

Sin^B deg day"*. 


Rotation rates of individual sunspots 
were determined only from one-day differences 
of longitude mesurenients. In order to identi¬ 
fy individual sunspots from' one day to the 
next using a computer program, which is a 
difficult task even when done by eye on a 
small sample of spots, we first grouped the 
spots into sunspot groups by proximity. 
Then we determined the identification of 
groups from one day to the next. The final 
step was to match spots within the group 
on two consecutive days using a technique 
which maximized the number of spot coincidences 
within a number of possible relative orienta¬ 
tions (excluding rotation of the group) in 
a two-dimensional grid. This method turned 
out to be insensitive to the assumed rotation 
used to.pick the group associations. 

Altogether nearly 36,000 groups were 
identified on the next day, and from these 
96,283 sunspots were identified as next-day 
returns. The method is quite conservative 
as i.s shown by- the fact that only 29% of 
spots were defined as returning the next 
day, and 10.5% of the groups seen the following 
day, contained no spots defined as returns. 






2.G5 

\ 



/ 

2.70 

\ 

\\ 


MI Ml SON 

twm.ii- 

4 

2 75 

V 

\ 

/ 

7; _ 

' 

200 

> 'A 

f 

1 



V 

CHOUMS - -1 1 


2.05 


Va 

' 

F--J -.*U SPOTS 

2.90 


, i% V 



2 95 

-30 

•20 -tO 

O 10 20 

JO 




LATITUDE 




Fig.1. Rotation rates ior all sunspots and for 
sunspot groups are plotted for 5-degree 
latitude zones. These results rofor to 
the Interval 1921 to 1902. Also plotted 
is the Doppler sunspot rate deter ml nod 
from the Mount Wilson data .for the 
interval 1967 to 1982 (Howard et a I., 
1 903), corrected for a small dispersion 
error (Snodgrass at al,, 1984), The Dopplqr 
dala are plotted from the smooth series 
expansion. 


Figure 1 shows the rotation rate 
of sunspots and groups as a function of latitude 
Also shown for comparison is the Doppler 
rate determined from the Mount Wilson magneto¬ 
graph data for the interval 1967-1982 (Howard 
eta 1.1983). 






The Rotation of tho Sun from Mt. Wilson Sunspot MuasurBmenta 
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For the largest sunspots (area > 
15 millionths of the hemisphere) the rotation 
may be expressed 

w ■ 14.282 (±0.018) - 2.636 (±0.169) X 

Sin 2 B deg day" 1 . 

For medium sized spots (5 < area < 15) 

to = 14.439 (±0.010) - 2.616 (±0.088) X 

Sin 2 B deg day" 1 , 

and for the smallest spots (area < 5) 

to = 14.549 (±0.005) - 2.869 (±0.051) X 

Sin 2 B deg day". 

The rotation rates of the sunspots 
of various sizes as a function of spot size 
are shown in Figure 2. 



Fig.Z. Tho rotation rate of spots of various 
sizes as a function of latitude. The 
zones are 5 degrees In width. 

Variations of the average sunspot 
rotation rate are seen which correlate with 
phase in the activity cycle. There is strong 
peak in rotation near sunspot minimum and 
a weaker peak near sunspot maximum. A similar 
pattern is seen in Mount Wilson Doppler rota¬ 
tional velocity data, which is available 
starting in 1967. There are at times small 
phase shifts between the peaks seen in the 
two data sets. 


The amplitude of the sunspot rotation 
variations is about 0.5$ day" 1 or about 3- 
4$ of the average equatorial rotation rate. 
The amplitude of the Doppler variations is 
nearly double that values 0.6 - 0.7 deg day 1 
or 4-5$ of the average equatorial rate. 


6. Discussion 

The variation of rotation rate with 
spot size is a puzzling result. Ward (1966) 
found a similar result for spot groups of 
various sizes using the Greenwich data. 
The difference between the largest spots 
and the smallest spots in this study is nearly 
2 %. One possible explanation for this effect 
is that all spots rotate at the same rate 
at some depth where they are linked to sub¬ 
surface flux ropes, but that the larger spots 
experience more viscous drag than do the 
smaller spots due to their larger cross-section 

as they move through the photospheric plasma. 
We tested this hypothesis by examining rates 
during years when the spot-plasma rate differ¬ 
ence was large and years when this difference 
was small. We could find no significant 
difference between the rotation rate differ¬ 
ences of large and small spots in these two 
samples. This implies either that viscous 
drag is not a significant factor in determining 
the rotation rates of sunspots or that the 
viscous drag is the same for spots of all 
sizes. 

In recent years there has been some 
debate about whether or not the sunspots 
and the photospheric plasma rotate at the 
same rate. Although this analysis does not 
give new results concerning the Doppler rates, 
it does not make it evident that the plasma 
cannot rotate at the same rate as all the 
spots. 
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In abso lute terms, our rotation 
rate for spot groups Is about \% slower than 
that found by Word (1966) and coincides closer 
than 0,1% with that found by Newton and Nunn 
(1951). The reason for these agreements 
is that Ward's results refer to all spot 
groups from positions that were not area- 
weighted (Greenwich). These results, therefore 
tend to be influenced by the smaller spots, 
which, as we found, rotate faster than the 
larger spots. Newton and Nunn used single- 
spot groups, which tend to represent large 
spots. Since our group positions are area- 
weighted, they represent the rotation of 
the largest spots, hence the agreement with 
Newton and Nunn is good. 

The cycle-related amplitude variations 
of the sunspot rotation rate imply periodic 
angular momentum exchange with deeper layers 
in the sun. This exchange cannot take place 
with other latitudes at the solar surface 
because such large amplitude fluctuations 

are not seen in the Doppler data (Howard, 
1976). Nor can the momentum be exchanged 

with the gas above the photosphere, because 
the inertia of this material is so low. 

We presume that the magnetic flux tubes of 
the spots are linked to deeper layers which 
share more or less in the same phase the 
variations seen in the photosphere. Alternati¬ 
vely the deeper layers at which the spots 
are linked may show no variations, and the 
surface rate of the spots may be determined 
by the drag imposed by the surface plasma. 

In any case the rotation variations 

may be considered tg represent yet another 
mode of torsional oscillation of the sun. 
These motions represent a somewhat higher 
amplitude torsional oscillation than those 
discovered earlier (Howard and laBonte, 1980; 
LaBonte and Howard, 1982). The sun seems 

easi ly excited in torsional inodes. 
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DISCUSSION: 


SIVARAMAN: When the spot is near the limb 

there will be foreshortening. This could 
result in the over estimation of longitude 
by more than 1/2°. How did you take care 
of this? 

HOWARD; We did not measure spots that were 
more distant than 60° from the central meridian. 


SIVARAMAN: You need torsional oscillations 

to explain your spot rotation. How will 
this fit in with the rotation rate above 
the sunspot zone? 

HOWARD; From the Doppler measure we find 
that the high latitudes share in the rotational 
modulation of the sunspot zones. 

GOKHALE: Is the variation of the rotation 

of spot groups in agreement with torsional 
oscillation from Dopplergram? 

HOWARD; No. The effect described has an 
amplitude several times that of the torsional 
oscillation. We have made a preliminary 
search for the torsional oscillation effect 
in our data with no success. 
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BAGARE: Would it. bo possible t.o study from 

your data, the occurence of retardation in 
proper motion of the spot groups as reported 
by Waldmeir? 

HOWARD; Probably yes, we have to analyse 
the data in order to look for these. 

VINOD KRISHAN: Is the missing penumbrae 

of some of the spots due to the overlapping 
of different spots in a sunspot group? 

HOWARD: I don't believe that the reason 

for missing penumbrae is known for certain, 
It is most likely due to the magnetic configu¬ 
ration within the spot group, which may or 
may not involve a nearby sunspot, 

BALASUBRAMANIAM: Is there any correlation 

between the area of the spot and the latitude? 
How do you account for areas of merging spots? 
Once the merger/coalition starts, is there 
any relation between the merger time and 
rotation rate? 


HOWARD: We have not looked for such an effect 
yet. When comparing the position of individual 
spots from one day to the next, we allow 

only one agreement for a spot. Thus, when 
we try to lineup second of the merged spots 

on the next day, there is no spot to associate 

It with. We have not studied this aspect 

as yet, 

ALEEMj Could you please explain the modern 
method of measuring the area of sunspot? 

HOWARD: Our techniques is designed for speed 
of measurement rather than accuracy. We 
place the cross-hair of our measuring instru¬ 
ment at two positions, which are digitized 
and stored for later analysis.' These two 
positions are chosen so that 1) one arm of 
the cross-hair points to the disk center 
and 2) the area of the quadrilateral defined 
by the two cross-hair positions approximate 
that the sunspot umbra. We estimate our 
accuracy as about 30$ for a single 1 umbra! 
measurement. 
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LINE ASYMMETRY AND MAGNETIC FIELDS 


Wi 111am Livingston * 
National Solar Observatory 
Tucson, Arizona 


1. Introduction 

The measurement of the asymmetries 
of Fraunhofer lines shows promise as a diagnos¬ 
tic of convection at the surface layers of 
the sun and stars. We discuss today the. 
observational evidence that line asymmetry, 
and by inference granular convection, is 
inhibited by surface magnetism. 

2. Origins of Line Asymmetry 

St’. John, working at the Mt.. Wilson 
150-ft tower in the first decades in the 
1900s, discovered what he termed a 'Limb 
Effect'; weak to moderate strength Fraunhofer 
lines display a red shift near the limb, 
and the amount of this shift is greater for 
the weaker lines. This 'Limb Effect* can 
now be understood as a consequence of granula¬ 
tion velocities and as such, is more correctly 
thought of as a "Centre-disk effect". High 
resolution spectra at di$k centre indicate 
that the bright granules are rising relative 
to the darker intergranular space, the material 
of which is falling. Low resolution spectra, 
which average over many granules, then possess 
a brightness weighted blue shift at disk 
centre which disappears at the limb because 
the motions are predominately radial. 

X ' 

Operated by Iho Association oi Universities 
for Research in Astronomy, Inc. under contract 
with the National Science foundation. 


A 

A solar granule may be considered as 
a convection cell. In the centre is a rising 
column of hot gas moving slowly from deep 
down, accelerating near the surface to conserve _ 
mass, radiating energy and thus cooling, 
then turning over and flowing down around 
the cell boundaries. A moderately strong 
absorption line, such as Fel 5303 A, is formed 
throughout the granular layers. As a result 
the extreme wings of the line, which sense 
the deepest slow moving layers are relatively 
undisplaced compared to its main body. Like¬ 
wise the line core is undisplaced because 
it is formed well above the granular cell 
where all outward motion has ceased. A plot 
of the bisector of such line has a resulting 
C- shape. Displacement of the extreme of 
the middle part of the C from its toe can 
be taken as a measure of convective vigour. . 

Distortion of the C - shape can arise 
from line blends, atomic fine structure and 
an unsyinmetric instrument profile (Dravins, 
et al. 1981). All Fraunhofer lines are in 
fact blended to some degree. By carefully 
selecting lines that are minimally blended 
and then averaging many comparable lines, 
together we obtain a meaningful bisector. 
Atomic fine structure, is avoided by confining 
ourselves, to iron lines. Finally the 1-meter 
Fourier Transform Spectrometer (FTS) has 
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a demonstrated symmetric profile and exhibits 
very low noise so far as line bisectors are 
concerned. 

The stochastic nature of granulation 
is dealt with by spatial averaging. A field 
of 2 arc-minutes is useful for defining the 
variation of line asymmetry across the solar 
disk. Ultimate averaging is the entire disk 
(10 5 -10^ granules). The internal scatter 
in full disk FTS bisectors corresponds to 
about 2 m.s"*. 


3, Influence of Magnetic Fields on C -shapes 

Observational evidence - the 'abnormal 
granulation' near sunspots - plus theoretical 
prediction . prompts us to investigate how 
C - shapes are affected by ambient magnetic 
fields. FTS observations at 2 arc-min resolu- 
ion made in plage areas and in quiet areas 
for comparison at the same limb distance 
indicate that the bisector amplitude (wing- 
core) is always reduced in the plage. The 
only exceptions to this rule occur very near 
the limb where the amplitude may or may not 
be reduced. 


4. Variabi lity of Line Asymmetry with the 
Solar Cycle 

Two data sets have been, examined for 
the temporal variability of the C- shapes 
(Livingston, 1982). The profile of Fel 5250,6 
A has been studied for the years 1976-1981 
using material acquired with the 13.5-m spectro¬ 
graph. Because such spectrometer observations 
are relatively noisy ( ~ 50 m.s" 1 ) many days 
must be averaged together, but when this 
is done the position of mid- C relative to 
the core is seen to have diminished by r ?m 
A over these years of rising activity. Beginn¬ 
ing in May 1980, FTS profiles became available. 
The noise in £ - shapes from the FTS is much 


less ~ 5 m.s' 1 . FTS archives for several 

strong Fe lines indicate a continuation of 
the trend found with the spectrograph through 
1982. •■■During 1983 the C - shapes began to 
strengthen. 


5. Comment on the Pic-du-midi Results 

Macris and Kosch (1983) examined the 1 
very best white light granulation photographs 
obtained in recent years at Pic-du-midi. ; y 
The average size of the granulation was deduced • 
by measuring the average intergranule distance. 
Tentatively they find that granule size becomes 
smaller with increasing solar activity. 

6. Conclusions 

Full disk line asymmetry which is a . 
signature of granular convection has been 
observed to be diminishing with'time until, 
recently when the trend apparently reversed. 

A comparison of low spatial resolution' plage 
vs quiet regions- on the solar disk suggests , 
that surface magnetic fields may inhibit 
line asymmetry. The average size of granula¬ 
tion likewise is reduced by solar activity 
according to French observers. We conclude 
that magnetic fields, expanding over the 
sun's surface from active regions, may interact 
globally with the solar granulation to signifi¬ 
cantly alter the transport of energy through 
this outermost layer. ■ y 
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DISCUSSION! 

BHATTACHARYYA: What is the data from which 

the Trench have concluded that the granule 
size depends on phase in the activity cycle? 

LIVINGSTON: They have studied a few (perhaps 
half a dozen) superior white light photographs 
taken at Pic-du-midi. Their method involves 
finding the average distance between granules 
near disk center in a quiet region. 

BALASUBRAMANIAM: Are there any changes in 

granule size close to spot penumbra? 

LIVINGSTON: Although the reference escapes 

me, I am certain that in the vicinity of 
active regions, and therefore near penumbrae, 
one finds "abnormal granulation" which is 
fragmented compared to normal granulation. 


VENKATAKRISHNANi Can this abnormal granulation 
be an extension of the Inhibition effect 
in C-shapes due to magnetic fields? 

LIVINGSTON: Yes, I assume it is the same 
mechanism. However, to account for the full 
disk changes we observe the effect most to. 
be global. 

HOWARD: Do the changes in asymmetry seen 

in integrated light come from the individual 
plages on the disk? 

LIVINGSTON: A good question which I am not 
prepared to answer quantitatively. As mention^ 
ed above, I have been thinking more of a 
global, rather than active region Influence. 
We do not know that the observed changes 
in line strength, equivalent width, for example, 
cannot be caused juSt by individual plages. 
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OBSERVATIONS OF THE EVERSHED EFFECT - PAST AND PRESENT 
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There exist good reviews on this topic 
by Schroter (1967) and Maltby (1964). The 
purpose of the present one on this occasion 
is to bring out fully the pioneering work 
of John Evershed in all its aspects in this 
field and also to update the earlier reviews. 

John Evershed (1909a,b,c) while analy¬ 
sing the sunspot spectra obtained by him 
at Kodaikanal for pressure determination 
in spots observed that the spectral lines 
show displacements in the penumbral region 

of sunspots. This was in fact a remarkable 
observation since spectroscopic observations 
of sunspots were started several years earlier, 
but line displacements were not detected 
previously in the spectra of the spot penumbrae. 
The first plate showing this line displacement 
was exposed on January 5 and 7 in 1909 - 

75 years ago - on a spot at latitude 9°N 
and longitude 31°W of the central meridian 
covering a region 4650 A to 4790 A with a 
dispersion of 1.08 mm A" 1 . He attributes 
the unusual amount of detail visible 

in the spectrum as due to the "clearness 
and steadiness of the air at the time of 
the exposure". . A preliminary measurement 
of ten of the best defined lines showed a 
mean displacement of 0.027 A (or 0.86 km 
sec - *) at the outer edges of the penumbra 
indicating a receding velocity on the limbward 
. side of the spot. He found an approaching 
velocity of the same magnitude on that part 
of the penumbra closer to the centre of the 

disc,. Continuing his investigations with 


150 spectra representing seven different 
spots in the northern hemisphere and four 
in the southern he established the following 
results: 

1. All spots whether in the southern or 
northern hemisphere show the same amount 
of line shifts when at the same distance 
from the centre of the disc and the sense 
of displacement is to the violet on the 
side of the penumbra towards the centre, 
of the disc and to the red for the limbward 
penumbra. 

2. These displacements are maximum for spots 
between 30° and 60° longitudes and dis¬ 
appear when they are within 10° of the 
centre. 


N 



Disc 


fig.1. Radial and tangential positions ot the 
spectrograph silt (or observing Evershed 
flow in sunspots. 
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3, He found maximum line displacements for 
a radial slit position and no displacement 
for a slit normal to this position (tan¬ 
gential slit, Figure 1). 

These observations convinced 
Evershed that the displacement of lines was 
solely due to the Doppler shift caused by 
the movement of the sunspot gases directed 
radially and horizontally outwards from the 
centre of the spot. This radial outflow 

of solar plasma from sunspots is now known 
as the s Evershed effect. The absence of the 
tangential velocity ruled out any rotational 
motion of the gases which conflicted with 
the concept of motion of gas within the spot 
by Hale (1908). Evershed (1910) continuing 
his studies using absorption lines of different 
strengths found out that strong lines (Rowland 
Intensity 5-6) show lower velocities 
( ~1.5 km/sec"*) while faint lines (Rowland 
Intensity 0-1) show higher velocities (2- 

3 km sec"*), He also established the monotonic 
Increase of radial velocity from the umbra 1 
boundary with a maximum near the outer boundary 
of the penumbra. . Although his earlier studies 
(1909a,b) showed that the motion abruptly 
ceases near the periphery of the penumbra, 
his subsequent observations (1916) gave evi¬ 
dence of the extension of the motion almost 
8000 km beyond the penumbra 1 boundary. From 
his observations using the Call H & K lines, 
he showed that there, is an inward flow of 

matter in the spot with a velocity of 
~1.8 km sec"*. Evershed also detected tangen¬ 
tial component of the velocity of ~3.6 km 
sec"* in five of the spots he studied and 
the direction of rotation was counter clockwise 
for some spots while clockwise for others. 

The voluminous work. of St. John (1913) 
using the 60-f.t solar tower at the Mt.Wilson 

i Observatory, of measurements of the Evershed 
velocities of 500 selected lines representing 


different elements and a large range of excit¬ 
ation potentials confirmed Evershed's earlier 
findings on faint metallic lines. For the 
faint lines he observed an outflow of materials 
from the spot into the undisturbed photosphere, 
whereas the strong lines (A 1 and Fe lines) 
showed no Evershed velocity. For still strong¬ 
er lines (Nal 0^, Dg, and Mg b^) and for 
chromospheric lines there was a reversal 
in the direction of motion indicating a radial 
flow of material into the spot (Figure 2). 



Fig.2. Schematic representation o( the height 
dependence ol tho Evershed How. Tha 
arrows towards tho right represent out- 
itow and the arrows towards the loll 
represent an Inflow of material. 

This would imply that the material depletion 
in the umbra due to outflow at the lower 
levels, is replenished by inflow of material 
at the chromospheric level. Evershed's (1910) 
attempts to detect motion in the umbra using 
the photospheric lines showed surprisingly 
a descending motion of~0,4 km sec"*. Abettl's 
(1932) observations during the sunspot maximum 
of 1926-1934 showed that the radial velocity 
varied from 0.2 to 6 km sec"* and in addition 
there were irregular tangential components 
'In the range 0-5 km sec"*. He also detected;: 
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that the maximum velocity occurred at the 
umbra-penumbra boundary and not at the outer 
edge of the penumbra, which contradicted 
Evershed's findings,. But Abetti's samples 
belonged to complex spot groups, whereas 
the measurements of Evershed and St. John 
pertained to single and rather inactive sunspots. 
This shows that there may be a clear distinct¬ 
ion in the pattern of the flow between the 
two varieties of spots. 

Michard (1951) was the first to relate 
the outflow velocity with the size of the 
spot and its magnetic field. He also establ¬ 
ished an empirical relation for the centre - 

to-limb variation of the velocity (V A ) 

r, o 

in terms of the maximum radial velocity 1 
(V„ „ ) and the angle (o) between the line 
joining the observer and the direction normal 
to the solar surface by the relation, 

V r,0 ■ V r,o <l-°.8 Sine) 

Later Schrtiter (1965a) obtained a similar ’ 
relation 

Ve * v r ,0 (1-0-3 Sine) 

This observed centre-to-limb variation confirms 
the height dependence of the velocity suggested 
by St.John (Figure 2), Kinman (1952,1953) 
in his classical work on this topic derived 
the spatial distribution of velocities over 
a sunspot from the spectra obtained at the 
Oxford solar telescope. Assuming cylindrical 
symmetry about the spot centre he derived 
the radial, tangential and vertical velocity 
components from the mean sight-line velocities 
at many locations within the spot. His Inves¬ 
tigation covering one large and three small 
spots showed that the maximum radial velocity 
attained ( Increased linearly with umbra 1 size. 


but the time to reach this maximum velocity 

was Independent of the size of the umbra, 

The interesting point is that the radial 

velocity component \l showed an Increase 
-1 " 

from ~ 1km sec at the edge of the umbra 
and reached a maximum value of~2 km sec" 1 
at the middle of the penumbra and finally 
reached zero in the photosphere well outside 
the spot. But according to Brekke and Maltby 
(1963) the radial velocities end abruptly 
at the outer boundary of the penumbra. 
Schroter's spectra confirm Maltby 1 s findings, 
but all these as we shall see soon depend 
to a very large extent on the spatial resolu¬ 
tion employed for the observations. The 
spot velocity fields derived by Holmes (1961) 
for two spots - one large and another small 
- using the non Zeeman line Fe I 5576.1 A 
showed a maximum radial velocity at the edge 
of the penumbra with an extended flow reaching 
zero well beyond, in the photosphere. But 
the ' maximum velocity for the large spot was 
only ~ 0.6 km sec" 1 as against the velocity 
\ of 3.5 km sec" 1 expected from Kinman.-s linear 
relation connecting maximum velocity and 
spotsize. Also for the large spot, the velo¬ 
city-reaches out far into the photosphere 
before it becomes indistinguishable, From 
the measurements of the radial, tangential 
and vertical velocity components for a spot 
at six disc positions, Servajean (1961) showed 
that the maximum radial velocity in spots 
increased with increments in y ( = Cos 0 ) 
values. His'interpretation Is that the higher 
layers have lesser radial velocity in spots 
as compared to the deeper layers. Holmes 
(1963) also detected a similar variation 
in the IL, as the spot approached the limb, 
but interpreted 1 this as the effect of scattered 
light in the telescope which becomes important 
as the spot is foreshortened at the limb. 
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High Resolution Observations 

Among these mention must be made of 
observations by Holmes (1961, 1963), Bhatnagar 
(1964) and Schroter (1965a,b) who used disper¬ 
sions in the range 5.0 to 8,0 mm/A for the 
study of the Evershed phenomena. 

Bhatnagar (1964) measured the sight¬ 
line velocities on a spot during its successive 
passage across the disc for various positions 
on the disc (ranging from p = 0.35 - 0.95) 
using non-Zeeman lines (Nil 4912 A, Pel 5576 A 
and Fe 1 5691 A). He used the Kodaikanal 120 
foot focal length solar tower telescope which 
has an image scale of 5.5 arc sec mm" 1 and 
the associated. 60 foot spectrograph at disper¬ 
sions of ~ 8 mm A" 1 . The radial velocities 
derived by him from the line-of-sight veloci¬ 
ties assuming cylindrical symmetry, steeply 
rose from the umbra 1 border reaching a maximum 
value of 1,5 - 2.0 km sec" 1 near mid penumbra 
and then declined gradually to almost zero 
far out in the photosphere. Ni 1 4912 A and 

Fe 1 5691 A showed similar velocities whereas 

* 

Fe I 5576 A had systematically lower velocities 
by about 30-35%. The velocities in all the 
three line progressively decreased from the 
centre to the limb, a fact that was noticed 
by Mi chard (1951) and Servajean (1961) earlier. 
This is, more pronounced for observations 
made close to the limb (Holmes, 1963). The 
straight line fit to the plot of the velocities 
for different ^positions show that the slopes 
for the two weak lines 4912 A and 5691 A are 
nearly the same, while the slope is flatter 
for the strong line 5576 A. This makes one 
fee 1 that the velocity dependence with height 
in the atmosphere may be real, although the 
scattered light would contaminate the estima¬ 
tion of velocities to a certain extent. 
Another interesting point in Bhatnagar's 
velocity measure is that they pertain to 


the successive passages of the same spot 
and the fact emerges that the radial velocities 
are higher in the fully developed phase than 
in the early phase of growth of the spot.. 

The existence and the numerical value 
of the tangential velocities are debatable 
points. There are as many observations in 
the literature showing the presence of the 
tangential velocities as there are pointing 
to their absence. Evershed (1909c, 1910) 

detected tangential components of 0.25 to 
0.35 km sec However in his later measure¬ 
ments (1916) he found no evidence for a regular 
tangential component, although irregular 
velocities were present at times. Kinman 
(1952) believed that the tangential component 
was random in nature and in general was within 
the errors of measurement although at some 
locations within the spot the value was as 
much as -0.5 t 0.2 km sec" 1 . But his 
neglect of the tangential component in the 
interpretation of his later measurements 
(1953). appears prejudicial , Abetti 1 s (1932) 
measurements on 26 spots show irregular tangen¬ 
tial component with values varying from 
Q to 5 km sec 1 . Measurements of Servajean 
(1961) and Bhatnagar (1964) are indecisive 
and the errors in the latter's measurements 
are rather large.making a realistic conclusion 
on the tangential components not possible 
Systematic and precise measurements with 
high spatial and spectral resolution are 
called for to settle this question. 

The use of high spatial resolution 
and spectral resolution hot only made the 
velocity measurements more precise and error 
free, but also brought out the microscopic 
structure of the Evershed flow as motions 
in the penumbra! filaments as we!) as the 
phenomenon of asymmetry of spectral lines 
in the penumbra 1 regions. Bumba (1960) using 
the high dispersion spectrograph of the Crimean 
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Observatory detected the strong asymmetry 
in the line profiles within the penumbra 
which accompanies the Evershed flow. Servajean 
(1961) studied the changes In the line asymme¬ 
tries for different spot positions on the 
solar disc. Bumba termed this phenomenon 
the 'line flag' whereas Holmes (1963) who 
also studied the asymmetry called it the 
•line flare'. It is interesting to note 
that Evershed (1916) detected in his low 
dispersion spectra that the spectral lines 
develop a 'diffuse wing' in the direction 
of the Evershed displacement near the outer 
boundary of the penumbra. Similar diffuse 
asymmetric widening was also detected in 
the spot spectra at the McMath-HuIbert Observa¬ 
tory (1956). But its real significance was 
brought out only by the work of Bumba and 
Servajean. Under very favourable seeing 
conditions, Bumba was even able to resolve 
the broadening into a satellite line but 
Servajean, Holmes or Bhatnagar do not find 
any Indications of such a satellite line 
in their spectra. Bhatnagar (1964) in his 
photometric study has investigated the depen¬ 
dence of the asymmetry on the strength of 

the line, on the locations of the slit on 

the spot and on the spot's position on the 

disc. The asymmetry grows In magnitude from 
the edge of the penumbra and reaches maximum 
values around the middle of the penumbra. 
The asymmetry has the same sense as the Ever¬ 
shed flow for the three lines he 
used (4912 A, 5576 A and 5691 A) and probably 
therefore is independent of' the strength 
of the line. Also, the maximum asymmetry 
is a function of the spot's position on the 
disc, such that the spot near the disc centre 
has an asymmetry maximum smaller than for 
spots near the limb. Again stronger lines 
show less of asymmetry than the weaker lines. 


Bumba finds that the flagging in l)j_ and 
lines of Na I is opposite to the direction 
of the Evershed displacement and the direction 
of flag near the core of these lines Is oppo¬ 
site to that in the wings. Servajean (1961) 
has suggested that the asymmetry is a manifes¬ 
tation of the different velocities correspond¬ 
ing to the different strata in the line forming 
layers, whereas according to Bumba (1963) 
the upstreaming of material at the umbral- 
penumbral boundary which bends and becomes 
horizontal In the penumbra cause the asymmetry. 
The asymmetry can also result when the outward 
moving material in the bright fibrils of 
the penumbra cause a blue wavelength shift 
which is not balanced by the Inflowing matter 
in the dark fibrils of the penumbra (Schroter 
1965b). Such an asymmetry would give a C- 
shape to the bisector of the line profile. 
The C-shape and its change for different 
p positions have been investigated in great 
details by SchrSter (1965a) using the Freiburg 
spot spectra employing the equi densitometry 
technique. The potentialities of this techni¬ 
que in bringing out the finest details in 
the spectra are well illustrated in his paper. 
His (Schroter 1965a) figure 6 shows that 
the displacement Is small In the line core, 
and the asymmetry reaches a maximum value 
at a distance of 0.8 times the radius of 
the spot from the spot centra. The centre 
to limb variations of the Evershed velocities 
at a fixed location from the centre of the 
spot shows that the velocities increase towards 
the limb reaching a maximum at u a 0.4. He 
also studied the variation of the asymmetry 
using the line Fe I 5123 A within the penumbra. 
The core of the line shows small displacements 
which increases away from the core and again 
falls less steeply at the wings. According 
to him, the sunspot spectra obtained at high 
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spatial resolution show nofr'Orfiy' 5 line asymmetry 
but resolve into strong satellites' (as seen 
by Bumba), which appears diffuse at the photo¬ 
sphere - penumbra ' boundary and rather sharp 
on the other side of the line centre, 

.The Evershed flow in the photosphere 
-chromosphere interface region was investi¬ 
gated by Bones and Maltby (1978) from the 
velocities in the core and in the wings of 
the hfg line using high resolution spectra. 
The measured velocities ( ~ 1.8 km/sec) follow 
the pattern already known from earlier studies. 
Another front from which the motions in suns¬ 
pots were investigated was by using high 
spatial resolution photographs of sunspots. 
From the time sequence photographs of the 
Princeton stratoscope project Schroter (1962) 
was. able to detect the systematic motions 
of the bright knots in the bright penumbra 1 
filaments. These" knots which form part of 
the fi laments move outward from the centre 
of the spot with a mean velocity of ~ 2km 
sec" 1 . He also detected material motion 
in the opposite direction in the dark channels 
that lie in between the bright filaments. 
There could be little doubt that these proper 
motions ,of the bright knots represent the 
Evershed flow that we see in the spectra 
of spots. Some of Bhatnagar's spectra (1964) 
show continuum brightness fluctuations In 
the penumbral regions that represent the 
motions of the-bright knots, Maltby's (1975) 
results, on the velocities In the chromosphere 
derived, from H fl filtergrams obtained on 
either wings at AX + 0.25 A, + 0.5 A and 
+ 0.75 A with a resolution of 1/8 A provide 
a three dimensional picture of the Evershed 
flow and its variation with time. The flow 
which is concentrated in channels is directed 
towards the umbra. These velocity channels 
occur mostly in regions occupied by the 
dark fibrils, whereas the flow channels with 


velocity directed away from the spot appear 
bright as in the continuum (Schroter, 1962). 
One interesting point is that the apparent 
length of the flow channels differ at different 
heights - they ‘are long in the chromosphere 
(- 16 arc sec) and shorter at the lower levels 
(~ 7- arc sec). The flow channels which are 
in the form 1 of loops alter with a characteris¬ 
tic time of about 10 minutes which possibly 
also represents the characteristic time of 
perturbations of the fldw channels. 

Finally, I shall mention about the 
high spatial resolution two dimensional obser¬ 
vations of Sheeley and Bhatnagar (1971) and 
of Sheeley (1972) which have provided a final 
answer to some of the questions relating 
to the Evershed flow. For example the question 
where the Evershed flow terminates has been 
a subject of much dispute since the time 
of its discovery. They derived a series 
of Dopplergrams from a time lapse sequence 
of spectroheliograms obtained on either wings 
of the Zeeman insensitive line Fe 1 5434 A. 
These velocity field pictures showed clearly 
that the Evershed flow confines itself to 
the dark spoke like filaments that end in 

the rugged outer boundary of the penumbra 

and the magnitude of the Evershed velocity 
is -4-6 km sec" 1 . Besides this, they also 
noticed an extra penumbral horizontal flow 

with speeds of 0.5 - 1.0 km sec" 1 which exten¬ 
ded to 10,000 - 20,000 km beyond the penumbral 
boundary. Comparison of these Dopplergrams 
with the CN 3883 A spectroheliograms showed 
that this extra penumbral flow terminated 
at the facular network. This high spatial 
resolution study established that the extra 
penumbral flow is indistinguishable from 
the well known supergranular flow and is 

distinctly different from the Evershed flow. 1 

In conclusion 1 shall list the salient points; 
which require immediate attention and study.; 
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I, The relation between the Evershed veloci¬ 
ties and the size of the spot and the 
magnetic field configuration associated 
with It should be firmly established. 

II. We do not know how the motions are 
related to the evolutionary aspects 
of the spots? At what phase of the 
growth of the sunspots does the flow 
start? How are velocities related 
to the age of the spot? 

III. The problem of tangential velocities 
should be settled. 

IV. According to Galloway's model of the 
Evershed flow, the flow is in the form 
of pulses. It would be interesting 
to check for any periodicities or the 
pulse nature of the flow. 

V. How are the line asymmetries caused? 

VI. What is the nature of the velocity 
pattern in bipolar spots? 

Answers to these questions from studies at 
high spatial and spectral resolution should 
help us to understand fully the macroscopic 
structure of the Evershed phenomenon which 
even to this day remains as enigmatic as 
it was 75 years ago. 
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DISCUSSION: 

G0KHALE: Has anyone looked . for short time 
scale variations in Evershed velocity? 

SIVARAMAN: In the model of Galloway he shows 
that Evershed phenomenon is in the form of 
pulses* it would be worthwhile to obtain 
time lapse spectra and search for any temporal 
periodicities. 
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Abstract 

The role playod by convection in the formation of slnndoF magnetic flux tubes and 
in the dynamics of the gos within the tube is discussed in tire case of simplified models. 
Convection instability cannot drive systernntic downflows whereas convective buffetling 
of the tube con. The inclusion of hunt transport reduces the efficiency of convective 
collapse for the formation of strong fields. The implications of these *wo results for 
tiro solar magnetic flux lubes is pointed out. 


I. Introduction 

The interaction of velocity anti magnetic 
fields has been well studied in the context 
of magnetoconvection. The work of Parker 
(1963), Weiss (1966), Proctor and Galloway 
(1979), Galloway and Weiss (1981), all indicate 
that in a Soussinesq fluid the final result 
of the interaction of a circulatory velocity 
field and a uniform magnetic field is the 
formation of a structured magnetic field 
devoid of flow and a region of flow bereft 
of magnetic field. It is natural to ask 
the question whether such structures can 
be formed even in the case of compressible 
magnetoconvection. This question is particu¬ 
larly important in the light of the highly 
structured nature of the magnetic fields 
observed on the solar surface. 

One way of answering this question 
would be to assume a family of equilibrium 
structures and examine their stability as 
a function of some parameter characterising 


the fanil ly, e.g. the magnetic field. Stabi lity 
analysis of two dimensional structures (e.g. 
an axisymmetric tube with area of cross-section 
varying with height) is exceedingly difficult. 
The slender flux tube approximation, introduced 
conceptually by Defouw (1976) and established 
rigorously by Roberts and Webb (1978), helped 
in the linear stability analysis of such 
structures. Webb and Roberts (1978) showed 
that tubes with field strengths smaller than 
a critical value are prone to a convective 
instability. The critical value depends 

on the depth of the tube as well as on the 
superadiabatic temperature gradient. Spruit 

and Zweibel (1979) calculated the growth 

rates of the instability for tubes embedded 
in a realistic convection zone environment. 
Their critical field strength was » 1500 G. 
These values compare well with the observed 
field strengths. 

However, one does not know a priori 

whether the convective collapse of weaker 
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tubes would result in a hydrostatic tube 
or hydrodynamic tube within the time scales 
of interest in the solar context. Hasan 
(1983) calculated the nonlinear development 
of convective instability of a slender flux 
tube and obtained final hydrodynamic states 

independent of the initial value of the magne¬ 
tic field. In this paper we review some 

of our results of similar calculations for 
po'lytropic tubes, both for adiabatic (Venkata- 
krishnan, 1983) as well as non-adiabaic (Venka- 
takrishnan, 1984a) conditions. 

Furthermore, the magnetic structures 

observed in the photosphere are constantly 

buffetted by granulation and waves. We shall 
also consider the nonlinear response of magne¬ 
tic tubes to such perturbations. 

2. Nonlinear Development of Convective 
Instability with Slender Flux Tubes 

The various parameters that determine 

the evolution of a polytropically stratified 

slender flux tube are Po , the ratio of initial 
gas pressure of magnetic pressure, So the 
initial superdiabatic gradient and the bound¬ 
ary conditions. The basic equations for 
a slender flux tube were first rigorously 
formulated by Roberts and Webb (1978). These 
equations are: 

it (p/B) + (P ' ,/B) ■ 0, (Z.i) 

If + 4* + !i p + 3 = 0 (2-2) 

(L + vL] p .IE(L+ v l_\ p = 

\9t 9z ;p p 'at 9z p 

■ -(y-1) (V.F) (2,3) 

B 2 = 8n (p e - p), (2,4) 

where p, p , v , B are the pressure, density 


velocity and magnetic field inside the tube, 
is the pressure outside the tube and F 
is the energy flux. We shall first consider 
the case F = 0 corresponding to ■ adiabatic 

motions. 

2.1. Adiabatic Flow 

The equations (2.1) through (2.4) with 
F * 0 form a system of hyperbolic partial 

differential equations which can be cast 
in the characteristic form. We integrated 
these equations using a backward marching 

algorithm based on the method of characteris¬ 
tics (For details see Venkatakrishnan, 1984b). 
In this way we obtained the values of all 
the dynamical variables as a function of 
z for various instants of time. In this 

paper we shall mainly discuss the evolution 

of velocity and magnetic field of the tube 
at a few spatial locations. 

We applied two different sets of boundary 
conditions. In one set the Eulerian pressure 
was kept constant at both ends of the tube. 
These conditions could be more appropriate 
to closed magnetic regions of the sun (Venkata¬ 
krishnan, 1983). The tubes on which such 
conditions were applied will be called 'closed 1 
tubes In what follows. In the other set 
of boundary conditions the Eulerian pressure 
was kept constant at the base of the tube 
while the Lagrangian pressure was kept constant 
at the top of the tube. This type of boundary 
conditions would be more appropriate for 
open vertical tubes in the sun. These tubes 
will be denoted as ’open' tubes in what follows. 

The general behaviour' of velocity for 
the adiabatic case consisted of an initial 
slow increase till a time t and then a rapid 
increase beyond that time. The onset of 
the nonlinear phase was delayed for stronger 
initial magnetic fields (small values of 
3o)* It was apparent that t-*«> as $ 0 tended 
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to some finite value ®2.0 (Figure 4 of Venkata- 
krlshnan, 1963). Furthermore, the Intensifica¬ 
tion of the field was not seen t* be sustained 
during the nonlinear phase for. Bo 3 4.0, 5.0 
and 6.0. An .Interesting example of oversta¬ 
te Illy was seen for Bo 5 2.0 In the case of 
an open tube. When Bo = 4.0, however, it 
was seen that steady magnetohydrodynamic 
states are attained. For Bo = 6.0, the non¬ 
linear regime showed unsteady behaviour even 
unti 1 30 free fa 11 times. 

The following picture emerged out of 
these calculations. The magnetic field inhi¬ 
bits the convective collapse of flux tubes. 
For some value of the field, overstable oscill¬ 
ations set in while for weaker fields steady 
dynamic states are attained. Even this is 
not possible for the weakest fields studied 
here, and in these cases unsteady behaviour 
persists for large times. The actual values 
of field strengths demarcating these regimes 
depends on the boundary conditions, with 
open tubes exhibiting greater stability compar¬ 
ed to closed tubes. 


2.2. Effect of Radiative Heat Transport 

The energy equation for an optically 
thick slender flux tube is: 


dft ¥ lr5M“) <!t.+ v !7)p =(y-Dk 


8 


■at ■ 3z* 


3 

■at 
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4(Te-T) 


9 Z T 


r 2 
1 0 


9Z 


* * (fi > < 


1 9 . 9j3 9 , 

l 3T 3T 3p' 


InK], (2.! 


(For a derivation of this equation see Venkata- 
krishnan 1984a). The first term on the R.H.S. 
of the equation (2.5) represents the lateral 
heat exchange. The second term represents 
longitudinal heat diffusion for constant 
heat conductivity K, while the third term 
represents the influence of variable conduct¬ 
ivity. Here we shall discuss only the case 


of constant K, 

When lateral heat exchange alone was 
considered In an open tube, it. was seen that 
oscillations were produced. The amplitude 
of these oscillations Increased with decreasing' 
radius r Q of the tube. For r Q <0.5 the 

oscillations became so large that local 
pressure enhancements rendered the magnetic 
field imaginary and precluded further computa¬ 
tions. The amplitude of oscillation was 
seen to be drastically reduced for 3 0 = 4.0, 
as compared to the case for Bo = 6.0. When 
longitudinal heat transport was included, 
the frequency of the oscillations was seen 
to be doubled and the amplitude reduced. 

This perhaps could be due to the excitation 
of a new mode or harmonic. 

Thus heat transport was seen to shift 

the state of marginal stability to lower 
values of the magnetic field as compared 
to the adiabatic case. Wherever the instabili¬ 
ty was evidenced, it was only in the form 
of oscillations and not monotonic as seen 

in the case of adiabatic flow. 

3. Response of Tubes to External Turbulence 

Let us for a moment forget the considera¬ 
tions of stability and assume that a stable 
magnetostatic flux tube exists in the solar 
atmosphere. This tube cannot remain static 

since it is embedded in a highly chaotic 
environment. The .environment can influence 

the tube broadly in 3 ways. It can twist, 
bend or squeeze the tube. A twist cannot 
influence the dynamics of the gas inside 
the tube while bending and squeezing motions 
can. In what follows we report some results 
concerning the''response of this gas to bending 
motions of the UlSte (Hasan and Venkatakrishnan 
1980) as well the response to squeezing motions 
(Venkatakrishnan, 1984b). 
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3*1. Response of Squeezing Motions 

The equations (2.1) through (2.4) descri¬ 
bed the state of a tube with a time independent 
environment. If external pressure perturba¬ 
tions are applied then they would induce changes 
in the internal gas pressure and, therefore, 
create a flow of gas along the field. We 
applied two forms of external perturbations. 
In one form the fluctuations were oscillatory 
in time but decreasing exponentially with 
height. The response of the gas within the 
tube is maximum when the scale length of 
the perturbation is approximately equal to 
the distance travelled by a tube wave in 
one period of the oscillation (Figure 1). 



TIME 

Flg.1. Response of gas to squeezing motions 
on Lhe tube - effect of exponentially 
decaying fluctuations with different 
values of scale height H. 

Moreover, the resulting flow along the tube 
is oscillatory. When the scale length is 
smaller than the tube travel length, the 
velocity amplitude is smaller and the flow 
tends to drift into a down flow. 


The second form of the pressure perturba¬ 
tion was chosen to be oscl llatory in space 
as well as in time. In this case,the maximum 
response is seen to Occur when the. wavelength 
of the perturbation matches with the tube, 
travel length (Figure 2). This confirms 
the resonance predicted by Roberts (1979) 
on the basis of linear analysis. It must 



Fig.2. Same as In Figure 1 but for wavollke 
pressure fluctuations for various ratios 
X of wove to lube speed. 

be seen, however, whether the resonant build 
up of amplitudes leads to a shock wave or 
whether the resonance is destroyed at large 
times. 

For the solar atmosphere if one takes 
the dimensions of the scale height of granular 
velocity decay into account then it would 
seem from these results that for strong 
tubes the most preferred response would be 
an oscillatory flow.' 

3.2, Response to Bending Motions 

In order to study bending motions, 








23 


Convection ond the Phenomenon of KllogausB Magnetic Fields on the Sun 


it is convenient to consider the motion of 
Infinitesimal flux tubes with area of cross- 
section 6 A which is inversely proportional 
to the local field strength. For this case, 
we Imposed a lateral velocity for the field 
line and studied the resulting flow along 
the field. The equations relevant to this 
study were written in a frame of reference 
moving with the field line, These euations 
are (Hasan and Venkatakrishnan, 1980): 


dVs + v 
dt s ds 


" p ds - 9 cosS + 

+ v n v s f, (3.1) 


d.iUP + v d ...In p + dVs 
dt ds ds 



+ 



3Vn 
9n ’ 


(3.2) 


p/p r b constant, 


(3.3) 


where p, p , V g and G are the gas pressure, 
density, velocity along the field, and the 
angle made by the field to the vertical respec¬ 
tively. The coordinate 's' is the distance 
measured along the field and V n is the velocity 
normal to the field direction. Once again 
these equations were integrated by the method 
of characteristics. A'detailed description 
of the resulting motion for two assumed forms 
of V n is given in Hasan and Venkatakrishnan 
(1980). The main result is that gas is . seen 
to be accelerated in the direction of increasing 
V . Physically this effect is because of 
the centrifugal acceleration on the gas cons¬ 
trained to move in a curved path because 
of the non-uniform V n< 

This effect can be applied to the case 
of granule flux tube interaction. It is 
known that the rms granular velocity decreases 
with height (Durrant et. a 1.: 1979); Thus 
if’ the tube is bent /with ■ a corresponding 
lateral. ve locity deereas 1 ng. with height, 


then, a downflow would result. A numerical 
simulation of this effect was performed by 
Venkatakrishnan and Hasan (1981), The result¬ 
ing velocities (Figure 3) compared favourably 
with some observations (e.g. those of G.lovan- 
e 111 et. al, 1978). A detailed account of 
the limitations of such an approach of modelling 
the interaction can be found in Venkatakrishnan 
(1984b). 



Fig.3. Inter action ol granolas with magnetic 

flux tubes - evolution of spatial velocity 
profiles. 

4. Application to Solar Kilogauss Fields 

There are two ways in which one can 
consider these results in the context of 
the tiny kilogauss flux tubes found in the 
solar photosphere. If one considers the 
.dynamics of the gas within these flux tubes, 
then it would appear from ail the above results 
that convective collapse with heat transport 
and external pressure fluctuations both produce 
only oscillatory motion. It is only the 
buffetting of granules that can be expected 
to. drive the observed downflows in these 
tubes. The problem of gas supply exists, 

however, and there is a school of thought 
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which-even doubts the reality of the downflows.' 
However! the fact remains that the stochastic 
. jostling of flux tubes with a spatially statio¬ 
nary gradient in the rms amplitude would 
necessarily lead to a downflow. The question 
. only remains as to the exact magnitude of 
the downflow. Recent observations using 
the 1-m FTS at Kitt Peak National Observatory 
(Stenflo et al. 1904) do demonstrate the 
necessity of including velocity gradients 

in flux tube models. 

The second aspect concerns the problem 
of convective collapse. For polytropic tubes 
we have; shown that radiative heat diffusion 
shifts the state of marginal stability to 
lower values of the field and introduces 
oscillations of the field strength as well. 
It remains to be seen whether calculations 
for non-polytropic realistic tubes would 
yield , similar results. Such calculations 
would , put. better theoretical limits on the 
minimum field strength of the flux tubes 
and also provide the exact: amplitudes and 
periods of the oscillations, if any. 

References 

DcfoUwi R.3.1?76, Astrophys. 3. 209, 266. 

Durrani, C.J.Mattlg, . W., Nesis, A., Reiss, G. 

• ‘ and Schmidt, W., 1.9 7 9, Solar Phys. 61, 

' 251 • '■ . . 

Galloway, D. J,, and Weltis, • N .0., .1901, Astrophys. 
J, 243, 945. 

G.lovanelli, R.p., Livingston, W.C., and Harvey, 
.1,, 1970, Solar Phys. 59, 49. 

Hasan, S.S., 190 3, IA U Symp. No. 102, 73. 

Husun, S.S., and Venkotnkrfshnnn, P., 1980, Kodoi- 
kuna! Obs. Dull. Set.A., 3, 6. 

Parker, L.N., 1963, Astrophys. J. 138, 552. 

Proctor, M.R.E., and Galloway, D.J., 197V, J. 

ritiiti Mech. 90, 273. 

Robotic, D., 1979, Solar Phys. 61, 23. 


Roborts, B., and Webb, A.R>, 197 0, Solar, Phys. 
56 , 5. 

Spruit, H.C., and Zweibel, E.G., 1979, Solar Phys. 
62, 15. 


Stenflo, 3.0,, Harvey, J.W, Brault, 3.W., and Solankl, 


S., 1904, Astron. 

Astrophys. 131, 333. 

Venkatakr Ishnan, P., 1983, 

A, 135. 

J. A9trophys. 

. Astron. 

Venkatakrfshnan, p,, 1904a, 
(In press). 

.1. Astrophys 

. Astron. 

Venkatakrlshnan, P., 1904b, 
University. 

Ph.D. Thesis, 

Bangal dto 

Venkatakrlshnan, P., and 

Astrophys. Astron. 2i 

H asan, S.5., 
i 13 3. 

1901, J. 

Webb, A.R., and Roberts, 

B., 19.70, Solar Phys. 


59, 249. 


Woiss, N.O., 1966, Proc. Roy. Soc. London, Sei.A., 
293, 310. 


DISCUSSION: 

BALASUBRAMANIAM: What sort of observations 

can one make, to look for these oscillations 
in ; sunspots? How. do you think they will 
show up? 

VENKATAKRISHNANj Oscillations are already 
known to occur in sunspots. My work is aimed 
'at tiny flux tubes. Only observations of 
high spatial resolution can detect the oscilla¬ 
tions.of field strength. Velocity oscillations 
can be detected using the shift in the zero 
point of the V-Stokes profile .but one needs 
high spectral resolution, in this case. 

GQPALSWAMYt Can . your computations be used 
for some inferences regarding the three tier 
nature of solar convection zone? For example, 
in the middle convection tier a single strong 
flux tube gives rise to many thin flux tubes 
in the upper tier, where, the radiation becomes 
very important. Can the oscillation of the 
middle tier flux tubes be considered as the 
temporal shaking of the thin upper tier flux 
tubes? 

VENKATAKRISHNANj. One can always, think of' 
a scenario where . the hierarchy of magnetic 
structures is governed by. the hierarchy of 
convective cells. However this is beyond 
the scope of my work. Moreover my work is 
concerned with the dynamics of individual: 
tubes and can therefore say nothing regarding 
aggregates of tubes coupled to each other. 
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Abstract 

The origin o( solar magnetic field and the Mux tube nature of this field on the small¬ 
est observed scales are two key problems In solar magnetohydrodynamics. Despite the 
success of klnematlcal models based on the turbulent dynamo theory, continuation of 
thinking on alternative approaches Is Justified and may be even necessary. 

One such approach Is based on consideration of energy transfer In magnetogssdynaml- 

cal production of new magnetic fields near the base of the solar convection zone. It 

turns out that the field produced by a substantial transfer of energy from the flows will 

essentially be in the form of strong thin flux sheaths that may break Into ’cluster of 

flux tubes' similar to those observed on the sun. Quantitative similarity requires nonlinear 

4-1 3 

MHD Interactl,ons of flows with amplitudes HQ cm s on timescales HO s. Such Interact¬ 
ions may be provided by azimuthal mngnetoacouatlc waves at the longitudes where they 
cross one enothor, The required wave amplitudes can be provided by a small fraction 
of the outward onorgy Mux which must be blocked near the base of the convection zone. 
The resulting phenomenological model of the magnetic cycle provides for (I) two varieties 
of activity with distributions and phases similar to the XBP end sunspot cycle, (II) poleward 
migrations of largo scale fields, (ill) long lived coronal hole like magnetic structures during 
decline of the sunspot cycle and (Iv) relations between the decline of each cycle and the' 
amplitude of tho succeeding cycle. The basic postulate of this model remains to bo estab¬ 
lished from solutions of MHD equations. However if established, It will provide promising 
alternative to the turbulent dynamo model. 


1. Introduction 

Solar magnetohydrodynamics can be said 
to have begun with the discoveries of strong 
magnetic fields and Evershed flows in sunspots. 
The present scope of the solar MHD covers 
various problems concerning the following 
phenomena (e.g. Priest 1982); (1) sunspots, 
(11) quiescent prominences, (111) transient 
events like flares, active prominences and 
coronal transients, (Iv) coronal heating 


and coronal structures, (v) solar wind and 
interp lanetary magnetic field, (vi) solar 
dynamo and solar activity cycle, and (vi 1) 
the small scale magnetic flux tubes In the 
photosphere and chromosphere. The boundary 
conditions used in modelling of the first 
five phenomena are often based on an ad-hoc 
choice of photosphere (or some other surface 
like the Alfv^nic surface) as one of the 
boundaries. Secondly the often observed 
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Inhomogenelties and discontinuities In dlrec¬ 
ti ons transverse ' to the field are Ignored 
In most of the models of these phenomena. 

The progress In the modelling of these phenome¬ 
na will depend ultimately on the boundary 

conditions at the "innermost boundaries" 
of the magnetic structures (i.e. the depths 

at which the field lines originate) and also 
upon the physical conditions and processes 
in the individual flux tubes constituting 
the magnetic structures. 

Thus, the solar dynamo problem and 
the problem of thin magnetic flux tubes can 
be considered as two key problems in the 

solar magnetohydrodynamics. 

In this review, I shall briefly discuss 
the various approaches to the solar dynamo 
problem and present a somewhat detailed discu¬ 
ssion of a new phenomenological approach 
which will bring out a likely connection 
between dynamo problem and the problem of 
thin flux tubes. 

2. Behaviour of the Solar Magnetic Field on 
Different Length Scales and Time Scales 

Because of its key role in almost the 
whole of the solar MHO, a model of the solar 
dynamo should be ultimately consistent with 
the behaviour of solar magnetic field as 
observed (or observationally inferred) on 
different length scales and time scales. 
I have listed in Appendix I as many observed 
properties of the solar magnetic field as 
known from different observations which are 
relevant to the present discussion. The 
purpose of presenting this list is not to 
suggest any attempt to "explain" all these 
phenomena with a single theory, It is only 
to illustrate how vast is the amount of obser¬ 
vational information that a solar dynamo 
theory will have ultimately to deal with. 


3. Solar Dynamo Problem as the Solar 

Cycle Dynamo Problem 

Solar magnetic fields would not have 
posed any dynamo problem if the observed 
large-scale field could be understood as 
a manifestation of the remnant of a 'primordial, 
field' which must have been trapped in the 
sun during its formation., ' Theoretically, 
the presence of such a fossil field cannot 
be ruled out even though the sun must have 
gone through a fully convective state in 
the past (Piddington, 1982). However, such 
a fossil field, it present, seems to be too 
weak to be detected in the noise of the strong 
fields associated with the solar cycle who 
-se surface average near polar region does 
reverse every 11.y or so. Thus, the solar 
dynamo problem is essentially that of the 
dynamo which provides the fields that manifest 
as the solar cycle. 

4. Different Approaches to Theoretical 

Modelling of the Phenomenon 

A theoretical model of any phenomenon 
must ultimately satisfy the following primary 
and auxiliary criteria; 

C^s Consistency with the basic laws of physics 
(primary), 

(^; Consistency with the relevant well estab¬ 
lished observations (pYimary), and 
: Consistency with well established theories 
of other related phenomena (auxiliary). 

Depending upon the initial emphasis 
on or , three kinds of alternative appro¬ 
aches are possible, - 

Physical; (Initial emphasis on ) - Zero 
order model is constructed by solving 
as many equations representing relevant 
laws of physics as possible. 

Semi-phenomenological; (Mixed initial emphas¬ 
is) - Zero order model is constructed 
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by solving some simpler sub-set of basic 
equations and replacing other equations by 
plausible postulates which provide consistency 
with some Important observations. 

Phenomenological: (Initial emphasis on C^) 

- Zero order model is itself a minimal 
set of. plausible and mutually consistent 
postulates having "provisions for consis¬ 
tency" with as many observations as 

possible. 


5. Successes and Failures of Existing Solar 
Dynamo Models Based on Such Approaches 

The following solar cycle dynamo models 
have so far been attempted from the above 
approaches: 

Phenomenological: 

.Babcock's model. (1961), Piddington's 
(1971) model. 

Semi-phenomenologica1s 

i) Oscillator models: (e.g. Layzer et 
al. 1979). In these, the equations 

of continuity, momentum, magnetic 

non-divergence are solved whereas 
the energy equation is ignored, 

ii) Turbulent dynamo models: (e.g. Parker, 
1963 j St1x, 1976; Yoshimura, 1979} 

Krause and Radler, 1980). In these 
only the equations of magnetic induc¬ 
tion and non-divergence are solved 
whereas equations for continuity, 
momentum and energy are replaced 

by assumptions about the convective 
flows and turbulence. 

Physical: 

(Gilman and Miller, 1981). All essential 

MHD equations are solved, except that effects 
of turbulence are parametrized. 


Naturally, the physical models have 
the least a-priori 'provision for consistency' 
with the observed behaviour of the solar 
magnetic field and plasma motions on different 
scales. However, the importance of these 
models lies in showing how difficult it is 
to satisfy both the criteria, and C^, 

and how far away the present modelling atempts 
are from the ultimate goal of a satisfactory 
model. The relaxation of criterion Cj (ignor¬ 
ing the energy equation) in oscillator models 
has not yielded commensurate success in satis¬ 
fying the criterion C 2 (cf, Cowling, 1981). 
The turbulent dynamo models, which further 
ignore the momentum equation also, are relati¬ 
vely much more successful. Unfortunately, 
the turbulent dynamo models are open to several 
doubts and difficulties (e.g. Piddington, 
1981} Cowling, 1901) which are summarized 
in Appendix II. 

In the solar dynamo problem, the pheno¬ 
menological approach is most suitable since 
the abundance of observational information 
might be used to identify and Isolate most 
pertinent physical question's for further 
investigation. Thus, Babcock's model did 
provide for qualitative agreement with most 
of the then known behaviour of the solar 
magnetic fields on length scales $ R 0 and 
timescales ~lly (e.g. the polarity anu orient¬ 
ation laws of bipolar regions, migrations 
of Fields and butterfly diagrams). One great 
drawback of the model, however, was that 
the evolution of three dimensional topology 
of'Its magnetic fields - as required to account 
for the reversal of the global poloidal field 
was extremely unsatisfactory, In the framework 
of its own phenomenology, the magnetic buoyancy 
invoked to bring the toroidal flux ropes 
to the surface would more probably remove 
all the magnetic flux from the sun (e.g. 
Piddington, 1972), 
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An alternative phenomenological approach 
which I have developed since last several 
years (Gokhale, 1977, 1979a,b) Is intended 
to retain the advantages of Babcock's model 
and to eliminate the difficult situation 
that the magnetic buoyancy might eventually 
remove all the magnetic flux from the sun. 
This difficulty will be remedied if some 
MHD process keeps producing new magnetic 
flux in the sun on timescales £ 11 y, at 
an average rate needed to balance the loss 
due to magnetic buoyancy. 

6.1. The Timescale of the Most Efficienct 

MHD Generation of New Magnetic Flux 

Magnetohydrodynamical generation of 
a new magnetic flux implies transfer of energy 
from the plasma flow to the electromagnetic 
field at a rate exceeding the ohmic dissipation. 

For a flow of velocity y transverse 
to a magnetic field of intensity B q the time- 
scale *r of such an energy transfer must satisfy 

[p] Cu Z 3 M 4 i to] [u 2 ] [B 2 ] 

i.e. ■ 

M" 1 a ta] [B 2 ] [pf 1 , 

where square brackets denote the orders of 
magnitudes of the bracketed quantities, a 
is the electrical conductivity and p the 
density at the place of interaction. It 
is assumed here that there is no external 
electric field, so that the mean current 
density is ~ouB Q . The most efficient energy 
transfer would imply conversion of a substan 
-tial fraction of the flow energy and hence 
the corresponding timescale will be given 
by 


mV] w 4 

i .e. 

[T] 3 [o] 4 [B o r 2 [p] (1) 

6.2. Generation of Magnetic Flux Sheaths 
and the Resulting Clusters of Strong, 

Thin Flux Tubes 

For conditions near the base of the 
convection zone (o »10 8 emu, B Q k 3 G and 
p » KT 1 gm cm“ 3 ), equation (1) Implies? 

t s 10 3 s 

For velocities u « 10 4 - 10 5 cm s -1 (comparable 
to the velocity of solar rotation), the mean 
field intensity‘6* in the resulting electromag¬ 
netic structure will be w 10 4 - 10 5 G and 
the thickness of the structure will be as 
small as s 10 -10 cm. The other two dimensions 
of the structure will be comparable to the 
corresponding dimensions of the flow, which 
will be generally >> 6 . Hence, the newly 
generated structures will consist 
of electric current sheaths and associated 
magnetic flux sheaths. If the dimensions 

H, 1 to both u and 8^ Is « 10 10 cm (density 
scale height near the base of the convection 
zone), then the total flux in each magnetic 
flux sheath will be 

$ J3 B* H 10 21 -10 23 Mx 

Each such magnetic flux sheath may be expected 
to ‘crumple 1 sooner or later, into a ‘cluster 
of flux tubes 1 . The amount of magnetic flux 
4 > in an individual flux tube will be 

<t> ~ B* 6 2 s 10 18 -10 a Mx 

It is quite interesting to note that for 
u « 10 4 cm s' 1 the cluster will have a total 
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magnetic flux $ «10 Mx, which is comparable 
to the flux in a typical sunspot group, and 
each flux tube in the cluster will have a 
flux' tjj wlO Mx which is comparable to the 
flux in most of the individual flux tubes 
observed on the photosphere. 

6.3. What kind of MHD flows can attain 
velocities slO 4 cm s" 1 near the base of 
of the convection zone? 

From energy flux considerations alone, 
velocities ~10 4 cm s~* cannot be ruled out 
for convective flows near the base of the 
convection zone. However, the corresponding 
mass flux will imply too high velocities 
( ~ 10^ cm s" 1 ) near the surface, which are 
never observed. Hence, convective flows 
may not have velocities ~10 4 cm s“^ near 
the base of the convection zone, unless their 
mass circulation is somehow closed much below 
the observable layers, which is rather unlikely. 

Let us, therefore, consider oscillatory 

flows like those associated with 

waves. Out of the three types of waves 

with U XB q , the slow and Alfven modes are 

too slow to yield energy transfer on time- 

scales '-'10 s, On the other hand, the fast 

(l.e. magnetoacoustic modes with wavekngth 

~ R 0 do have timescales ~10 3 s. Among 

such modes the radially propagating modes 

cannot have amplitude *10^ cm s“^ near the 

base of the convection zone, since the 

corresponding radial energy ^ ux would be 

four orders of magnitudes larger 

than the total energy flux at the surface. 

Since the background field is poloidal, the 

only remaining kind of magnetoacoustic 

mode ll 1 are the azimuthal ones, and there 
— —o . 

seems no reason why these cannot attain 
amplitudes ^lO 4 cm s~*, 


6.4. Mechanism for Energizing the Magneto- 
acoustic Modes Near the Base of the 
Convection Zone 

Near the base of the convection zone, 
there will be a layer of finite thickness 
where the temperature gradient is superadlaba- 
tic but the superadlabaticity is below the 
critical value required for large-scale convec¬ 
tive instability. In such a layer, the minute 
fraction of the outward energy flux which 
cannot be carried further either by radiation 
or by convection will be "blocked". This blocked 
flux must go in energizing various waves 
including the azimuthal magnetoacoustic waves, 

6.5. Amount of Energy Required for Wave 
Amplification Upto 10 4 cm s“* 

For building up wave amplitudes s 10 4 
cm s" 1 the required average value of the 

blocked energy flux will be only' ~10® erg 

-1 

cm s l.e. only ss 0.1 - 1% of the total 
radial energy flux (also comparable to the 
mean value of the non-radiative flux deposited 
above the photosphere). 

6.6. Locations of Interactions 

Conservation of angular momentum requires 
that the azimuthal MHD modes be generated 
in pairs of oppositely propagating waves 
which keep on repeatedly crossing each other 
at a finite number of longitudes, where the 
steepest gradients of velocity and magnetic 
fields would develop. Therefore, it is in' 
the neighbourhood of such longitudes that 
the non-linear interactions satisfying equation 
• 41 ) may be localized. 

6.7. A Tentative Postulate 

On the basis of the foregoing discussion 
we arrive at the following tentative phenomeno¬ 
logical postulates 
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Near the base of the convection zone 
a minute fraction of the outgoing energy 
generates and amplifies oppositely propagating 
azimuthal magnetoacoustic waves and these 
waves produce strong, then, electric current 
sheaths and magnetic flux sheaths through 
nonlinear interactions near the meridians 
where they cross one another. The foregoing 
postulate remains to be verified by detailed 
magnetogasdynamical calculations. However, 
let 'us assume its validity for the present 
and proceed to see how one could develop 
from it a phenomenological model with provi¬ 
sions for agreement with several important 
observed characteristics of the solar cycle. 



Flg.l. A schematic diagram of the phanamono- 
logicnl model. Thick lines represent 
one ol the sovoral newly created magnetic 
"Mux rolls". Thin llneo represent the 
background megnotlc field "S". This 
background field (S) and the front portions 
(R) of the newly created "Flux rolls" 
emerge above the photoophore providing 
two varieties of activity with similar 
polorltieB but differing In phase, surface 
distribution and mognotlc flux per active 
region. Meanwhile the back portions 
of the flux rolls, represented by dashed 
linos, diffuse over all latitudes and 
longitudes to provide the reversed back¬ 
ground Held for the next cycle. 

For convenience of illustration 
only one newly creatod flux roll Is shown, 
and no attempt is mode to show effects 
of differential rotation. 


6-8. Consequences 

6.8.1. Reversal of the Global Field 

The newly produced magnetic structures 
must necessarily.be topologically closed. 
Thus the flux sheaths must be 'flux rolls' 
and the flux tubes must be 'flux loops'. 

Therefore, the reversal of the global field 
on an 11 y time scale can be achieved if: 
(i) the background field B Q and the similarly 
oriented front portions of the newly created 

flux loops rise upto and above the photosphere 
gradually on a timescale of 11 y and simultan¬ 
eously*, (ii) the back portions of these flux 
loops diffuse in the base layer to provide 
a background ofieId ~(-B q ) for the next cycle 
(Figure 1). 

6.8.2. The Two Varieties of Activity 

The global fiflld reversal in the above 
manner would imply that during any 11' y' cycle 
the activity at and above the photosphere 
will come from segment-by-segment emergence 
of the following two varieties of flux tubes 
(segments from both varieties having similar 
polarities): (1) The thin, frayed flux tubes 
of the background field which are’ present 

at all latitudes and all longitudes and whose 
current system and field lines' are confined 

to the outer parts of convection zone and 
the outer space, and (2) The front portions 
of the newly produced thick clusters of flux 
tubes which are present predominantly in 
the central latitudes and neighbourhoods 
of a few preferred longitudes, and whose 
current system and field lines are closed 
near the base of the convection zone. The 
first variety of segments will, in general, 
emerge earlier than the latter. 

Thus, the model of the field reversa 
automatically provides two varieties of solar 
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activity with the surface distributions and 
relative phases similar to the x-ray bright 
point and sunspot related solar activity. 
Schlissler (1980) presented a model based 
on a similar modulation of the thickness 
of flux tubes responsible for the activity 
without entering into the question of how 
the magnetic flux tubes are produced. However, 
in his model thicker flux tubes rise earlier 
which seems to contradict the fact that each 
new cycle starts with appearance of the small¬ 
est spot groups (e.g. Gokhale and Sivaraman, 
1982); 

6.8.3. The Background Magnetic Flux and the 
Activity Related Magnetic Flux 

Since the magnetoacoustic waves of 

amplitudes ~ 10 4 cm s" 1 are expected to produce 

21 

magnetic flux rolls of ~ 10' Mx each, the 
maintenance of the order of magnitude 
of the background magnetic flux requires 
that such a production occurs only once during 
each ~ 11 y or so, (e.g. in association with 
some 1 catastrophic transformation of the wave) 
if the azimuthal wave number of the wave 
is < 10. In such a case ('catastrophic 

model') each flux roll will have to emerge 
in ~10 - 10 segments in order to account 

for the ratio of that order between the acti¬ 
vity related flux and the background flux. 
If the required amplification is provided 
by the differential rotation as in the Babcock's 
model, then it will also account for the 
polarity orientations as in that model. 

Alternatively, if the magnetic flux 
rolls are produced during each wave period 
('persistent ■ .'production model 1 ) then the 
amount in the flux rolls and the rate of 
their emergence must vary during the 11 y 
cyc le in such a way that the order ' of rfiagni- 

tude of the background flux, is maintained 

and the number of segments give the 'right 


amount of the activity related flux. 

These details are yet to be worked 

out. 

6.8.4. Migrations of Large-scale Fields 

The emergence of the 'front' portions 
and the diffusion of the 'back 1 portions 
of the newly created fields will lead to 
poleward migration of the following polarities 
in the photospheric field along with the 
neutral lines separating these polarities 
from the background field (e.g. Howard, 1974* 
Makarov, et al. 1982; Makarov and Sivaraman, 
1983). 

6.8.5. Variations of the Sizes and Locations 
of the Long Lived Coronal Holes and Relations 
Between the Decline of One Solar Cycle and the 
Amplitude,of the Next Cycle 

During the declining years of the sunspot 
cycle’ the photospheric intersections of the 
new magnetic structures will provide deep 
rooted, monopolar, open field structures 
at middle ' latitudes which might account 
for the occurrence of long lived coronal 
holes at such latitudes during the declining 
phase of the solar cycle. • . 

Ultimately, towards the end of the 
solar cycle, these structures will yield 
monopolar regions of the reversed polarity 
near the poles. 

6.8.6. The Values of the Period 

The period of the solar cycle in such 
a model will be ultimately determined by 
the timescale of modulation of the rate of 
energy trapping near the base of the convection 
zone. As an .illustration it has been shown 
earlier (Gokhale, 1977) that a value of 
~11 y can be obtained (for a 'catastrophic 
model') by . assuming, a turbulent viscosity 
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6.8.7. Modulation of the Amplitude of the 
Cyc le 

Under 'normal 1 circumstances a relatively 
steady cycle may bo maintained by an approxi¬ 
mate balance among the rates of the following 
processes: (1) production, of magnetic flux 
near the b^e of the convection• zone, . (if) 
the rise of flu* tubes across the convection 
zone, (iii) migration and diffusion of magnetic 
flux across the Surface and (iv) the removal 
of magnetic flu* by magnetic buoyancy, reconne¬ 
ctions and solar wind. Excursions to and 
from 'abnormal' conditions like Maunder minima 
might be caused by some non-linear couplings 
between these processes, 

7. Conclusions 

I may conclude by saying that the theory 
of the solar dynamo Is far from the final 
solution of the problem, although the seml- 
phenumonologica 1 turbulent dynamo models 
and tho fu 1 ly phenomenalogica 1 new approach 
describes here continue to be the most promi¬ 
sing approaches. 

-The new phenomenology based on the 
postulated production of thin electromagnetic 
structures by locally non-linear interactions 
of magnetogasdynamical modes provides for 
agreements with many observed features of 
the solar cycle on timescales 422 y. However, 
it cannot be' considered • as / an .established 
model till the basic.’-.postulate .is verified 
.by mathematical /solution of - the - 'relevant 
magnetogasdynamic'a 1 equations, and the .para¬ 
meters of : the 'solution .are shown to - vary 
on timescales > 22 y in a way consistent 

with the observationaily derived behaviour 
of the solar cycle on such timescales. The 
dilemma of catastrophic versus continuous 
production of new field also needs to be 
resolved. 


Appendix I 

OBSERVATIONAL CONSTRAINTS 

Magnetic Cycle (Timescales ~ 22 y): 

1. Reversal of poloidal ( ~10 23 Mx) and toroi¬ 
dal fields after about every ~11 y or 
so. 

2. Phase of the reversal: ~1-2 y after sunspot 
maximum, 

3. Poleward migrations, of 'f polarities 
and magnetic neutral lines (filament 
channels). 

Sunspot Activity Cycle (Timescales <. 11 y): 

4. Periodicity in sunspot abundance with 
a period ~ 11 y (Historically the first 
'solar cycle phenomenon' discovered). 

5. Total flux: ~ 10 24 - 10 25 Mx. 

6. Distribution: 

aj central latitudes (* 40°) 
bj preferred longitudes (?) 
c) butterfly diagrams 

7. . Properties of the 'unit' of activity 

(viz. a spotgroup); 
aj polarity laws 

b) orientation laws,, 

c) typical flux ~l(r i Mx 


EAR and XBP Activity Cycle (Timescale & 11 y): 

B. Periodicity in Ephimeral Active Regions 
(EARs) and X-ray Bright Points (XBPs) 
with periods ~ lly. 

9. Phase: 

..a) similar to sunspot cycle for EARs 
b) opposite for XBPs 

10. .Total Flux: 

10 25 - I0 26 y.x 

.11/ Distribution: , ‘ . ■ ••' '• ; 

I. a) all latitudes : .(pole-to-pole) 

b ; ) ■ no' -..preferred -.. longitudes (except 
.the .tendancy to avoid those occupied, 
by-spot 'activity) ' ..' ,‘ 

12. Properties of the unit of activity (EAR/XBP) 

a) polarity' random but often -turning 
to normal for EARs 

. b) typical flux ~ 10 18 •• lO^Mx 

Magnetic field behaviour on smallest and 
medium scales length scales < R; 

13. 'flux tubes’ with B .~10 3 G; flux 
~lo¥: - 10 18 Mx. 


Solar Dynamo Problom 


33 


14. Clustering In active regions (upto 
~10 22 Mxi typical 10 21 Mx) 

15. 'Down flows' (?) total mass supply 

9 fi c 

10 gms per cycle. 

Coronal Behaviour on Timescale 

16. Variation of the shape and size of the 
corona with sunspot cycle. 

17. Variation in abundance, location and 
sizes of coronal holes (e,g, from geomag. 
data) 


Flows and Differential Rotation (length-scale 


10. Two varieties of differential rotation 
of magnetic structures. 

19. Variation between the two with the phase 
of the cycle. 

20. Torsional oscillations in the photospheric 
plasma rotation. 

Timescale < 22 y: 

21. Relations between (N) t,h & (N+l) th cycle 
(e.g. from geomag. data) 

22. Rise amplitude relations. 

Timescale £ 22 Y: 

23. Modulation on timescales ~ 00 y. 

24. Maunder minima and resurrection there 
from. 

Long Term Behaviour: 

25. Long time scale stability upto ~ 700 
mi 11 ioniums ago. 
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TURBULENT DYNAMO MODELS HAVE ACHIEVED AGREEMENT 
WITH MANY OBSERVATIONS. THEN WHY IS THE 
DYNAMO PROBLEM STILL UNSOLVED? 

1. Agreements with observations are not 
from a single set of input parameters, 

2. First order smopthening approximation 
requires |B 1 2 | < |B|. But observed thin 
strong flux tubes imply |B' [ « |B | if 
the tubes are densely packed (in which 
case w Lorentz-force's cannot be neglected) 
and |B I > |B' j otherwise. Here B is the 

.ensemble average and B' is the' random 
part of the field. 


3. Another basic assumption \>> <v> t Is not 
valid for the observed turbulence (Here 
X is correlation length, <v> Is r.tn.s. 
velocity of turbulence, t Is correlation 
time), 

4. No way of removing the small scale fields 
produced at the boundaries of the eddies, 

5. Requires dw/dt < 0. True near top of 
convection zone. But buoyancy removes 
the field completely (if at the base 
of convection zone where d u>/d t > 0 for 
stability, then ct effect is < 0). (Here 
w is angular velocity of rotation). 

6. When the dynamics is introduced, the 
field either grows indefinitely or contra¬ 
dicts migrations (butterfly diagrams) 
and the differential rotation (Gilman 
and Mi 1ler, 1981). 


TURBULENT DYNAMO THEORY CANNOT BE SAID 
TO HAVE SOLVED THE DYNAMO PROBLEM. CONTI¬ 
NUATION OF THINKING ON ALTERNATIVE APPRO¬ 
ACHES IS JUSTIFIED AND MAY EVEN BE NECESS¬ 
ARY . 
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DISCUSSION: 

GOPALSWAMY: Knowing the set of observations 

and the parameters to be introduced, can 
one do an optimization? 

GQKHALE: So far even the set of parameters 

has not been uniquely identified, and hence 
the question of the optimization of values 
does not arise. Turbulent dynamo experts 
must be working towards the final choice. 
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VECTOR MAGNETIC FIELDS IN PROMINENCES: Observations and Analysis 
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Abstract 

Measurements of the polarization of the ' Ho I multiplel in a prominenoe allow 

the magnetic fields associated with the prominence to be deduced. The amount and origin 

of the ilnear polarization are governed by the Hanlo offecl, the detailed; theory of which 
is applied to the Stokes components of the He I 0^ emission to derive the associated veotor 

magnetic fields. The dependence of the observed polarization on the magnetic field Is 

predicted from solutions of the statistical equilibrium equations for the He I D ^ mulllplet. 
Magnetic field vector solutions for each pair of llnoar polarization Stakes profiles correspon¬ 
ding to on obsorvatlannl point In a promlnonce are not uniquely determined. Nevertheless, 
physically correct solutions can usually bo identified with a high confidence level. .'■■.Fluid 
vectors are mostly close to horizontal end normal to the prominence surface, magnitudes 
ranging from 6 G to 60 G, with a mean of 24 G. 


1. Introduction 

The key to understanding the structure 
and origin of prominences lies in establishing 
their vector magnetic field configuration. 
Zeeman-effect magnetograph measurements allow 
only the longitudinal (line-pf-sight) component 
of the field to be measured with reasonable 
accuracy (±1 G) and with a moderate sensitivity 
down to ~1 G, while measurements of the trans¬ 
verse component have an estimated sensitivity 
of only 50 G (Leroy, 1979), a value apparently 

* Operated by the Association of Universities 
for Rosearch In Astronomy, Inc. under contract 
AST 70—17292 with tha National Science Foundation. 
** The National Center far Atmospheric Research 
is sponsored by the National Science Foundation. 


much larger than typical field strengths 
within quiescent prominences. Moreover, .It 
is now generally recognized that the validity 
of such measurements and the accuracy of. 
the corresponding derived fields, interpreted 
by the Zeeman effect can be severe jy limited 
by resonance polarization contributions. 
In the case of He I D 3 , for example, . thft 
observed circular polarization results both 
from the Zeeman effect and from level crossing 
interference in the Hanle effect (Land! Deg 1 1 
Innocenti, 1982). Also, radiative transfer 
effects in commonly used lines such as ^ 
substantially complicate the ..magnetic field 
interpretation. 
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Polarization due to resonance scattering, 
modified by the presence of a magnetic field- 
the Hanle effect - is a sensitive and accurate 
indicator of the vector magnetic field in 
a prominence, The Hanle effect occurs in 
a weak magnetic field when the Zeeman splitting 
is comparable to the natural width of the 
component sublevels* the resulting partial 
degeneracy of the level mixes the polarization 
properties of the Zeeman sub-levels. This 
mixing partially depolarizes radiation and 
reorients the plane of polarization of radiation 
undergoing resonance scattering. The extent 
of the depolarization depends upon the magnitude 
and direction of the field. From the amount 
of the depolarization and rotation, one may 
deduce the vector magnetic field associated 
with a prominence. 


He \ Q, 
Mulilplet 
587$ A 


-3r 


0.0 


A . .3 
MIA) 


.5 


fig-1. The relative in pgn I tu ties and spectral 
locations oi the transitions that form 
the He i D^ multipJet. 


The D 3 (5876 A) spectral profile has 
two well-defined components. The major (blue) 
component results from five transitions, 


3 3 D 2 -2 3 P 2 , 


and while the minor (red) component 


3V 3 >v 


3 3 D 1 -2 3 P 2 


3 3 

is a simple triplet, 3 D^-2 P Q , producing 
an extended wing in the profile. The peaks 
are separated by about 0.34 A. The relative 
magnitudes and spectral locations of the 
transitions are shown in Figure 1. The 
two components of the profile have different 
polarization sensitivities to magnetic field 
strengths (House and Smartt, 1979, 1982), 

and if .spectrally resolved and isolated by 
Gaussian curve fitting, two independent measure¬ 
ments of linear polarization (from the Q 
and U Stokes scans) are obtained, A histogram 
of linear polarization values, P(l), associated 
with the major component of the profile, 
representing 83 measurements of quiescent 
prominences, is shown in Figure 2a. Figure 



Fig.2. Histograms of the linear polarization, 
P, as observed In quiescent prominences 
corresponding to a) the major component 
of the profile (P(1)), and b) the minor 

component of the profile- (P(2)). 

2b is a corresponding histogram of the polari¬ 
zation, P(2), of the minor component re¬ 
vealing the significantly different distri¬ 
bution of P(l) and P(2), This difference 
in sensitivity produces some redundancy (for 
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perfect data) In the Information required 
to deduce the vector field. Further, Landi 
Deg 1 1 Innocentl (1982) has completed the 
theory of D 3 by computing the expected circular 
polarization profile .for a given configuration 
of the magnetic field vector. The observed 
circular polarization V profiles can then 
be interpreted. Since the resultant data 
from the V profiles are independent, amoigu- 
ities intrinsic to interpreting polarimetric 
data of prominences in D 3 can be resolved. 
However, ooservati ona 1 ly the V component 
is frequently noise-dominated and hence unin¬ 
terpretable, Additional advantages of 
are that it is both optically thin and relati¬ 
vely bright in prominences. Most Importantly 
however, the statistical equilibrium equations 
have been solved in sufficient detail (Bommier 
et al. 1981: Landi Deg 1 1 Innocenti, 1982) 
that D 3 may be Interpreted with confidence. 

Here we report on measurements of vector 

'<*• 

magnetic fields associated with a quiescent 


prominence, obtained by Interpeting spectrally 
resolved Stokes scans of 0 3 - For a more 
detailed description 1 ' 1 of this work, refer 
to House and Smartt (1982), Athay et al. 
(1983) and Querfeld et al. (1984). 

2. Observations ( 

The D 3 polarimetric observations of 
(the 'prominence discussed here were obtained 
on April 19-22, 1978, using the HAO/SPO 
Stokes-I polarimeter (Baur et al. 1980). 
A'*< set of Stokes profiles, coresponding to 
one observational point (3.5 arc sec x 10 
arc sec), represents 100 spectral scans recorded 
in about 20 minutes. The dispersion of the 
prominence in both H a and 0 3 on four consecutive 

days of observations, as well as on April 
23, is shown in the spectroheHograph images 
of Figure 3 (at approximately 14:00 UT), 
Note the significant evolution of the prominence 
over a four-day period. This prominence occurr¬ 
ed in the S-E, extending in longitude over 




April 19 20 21 22 23 


Fig.3. Spectrohellograph Images In and of a largo quiescent prominence of which polari¬ 

metric observations were obtained, April 19-22, 1 979. Approximately one hour after 
the last Image was recorded, the prominence erupted. 
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about 75 , and in latitude from approximately 
31° to 6 ? ♦ The long axis of the prominence 
was inclined by about 30° to the line of sight. 
The maximum height was ~ 10 5 km. This promi¬ 
nence erupted abruptly soon after 14:00 UT 
on Apri 1 23. 

A sample of a typical set of four D 3 
Stokes spectral profiles, Intensity, I , linear 
polarization, Q and U, and circular polariza¬ 
tion, V, is shown in Figure 4. The ordinate 



Fig.4. A typtosl set of Stokes profiles In D^. 

units are millionths of the central disk spectral 

brightness. The linear polarization degree 

(of the major component) Is about 2 % while 

the circular polarization is a factor of four 

smaller. Instrumental noise levels were below 
-3 

10 in polarization, Note that a two-slit 
magnetograph measurement of V would be very 
misleading. 

3. Data Processing 

Polar line trie observations such as these 
shown in Figure 4 are initially smoothed and 
fitted with Gaussians to obtain the relevant 
polarization parameters. As well, they must 
be augmented with supplementary information 


necessary for their interpretation. 

The polarization reference axis assumed 
in both the Bommier (1981) and Landi Deg 1 ' 
Innocenti (1982) statistical equilibrium 
solutions lies along the local radius so 
that Q and U must be rotated from the instru¬ 
ment coordinates, E-W geocentric, into the 
radial reference frame. In addition, the 
apparent height, h, of the observation above 
the limb and the dihedral angle, x? between 
the plane of the sky and the heliocentric 
radius through the observed point are necessary, 
h and x together give the true height, 
r = (R 0 + h) Cos x . The angle x is also needed 
to determine the disk-prominence-observer 
scattering angle. The measurement of x 
invariably uncertain. It i? best obtained 
by tracking a specific point on the prominence 
over a period of time. Such a sequence of 

observations, Ideally over several days, 
even if available, has limited precision 
due to the general evolution in form of all 
prominences even on a time scale of a few 
hours. 

All of the amplitudes for I, Q, U and 
V are normalized to the intensity to obtain 
blue and red values of Q/I, U/I, and V/I, 
proportional to the integrated strengths 
of the mu It1p let components. This results 
in four linear polarization values (Q/I) re j, 

(Q/Obiue- < U/I W < U/I >blue- wh1ch are 
used to obtain plausible, but usually ambiguous 

values for field strengths, |||, and direction 
(8 , 4 )), where 0 is the field polar angle measired 
from the local radius and 4 > is the azimuthal 
angle measured around the local radius, with 
<t>= 0° pointing west and 41 = 90° pointing north. 
The remaining data, (V/I) blue and (V7I) blue 
(where V is the derivative of V) are used 
to reduce the ambiguities inherent in the 
Q, U interpretation. 
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3.1. Hanle Effect Codes for He I D 3 

Bommier et a 1. (1981) and Landi Degl' 
Innocentl (1982) have Independently constructed 
Hanle effect codes for the magnetic field 
Interpretation of prominence Stokes profiles 
In He I Oy 

From the Stokes profiles, the observa¬ 
tional parameters of intensity and the corres¬ 
ponding two-dimensional linear polarization 
vectors for both the blue and red components 

of the mu 11iplet, as well as values for 

h and x > are employed in the interpretation. 

Essentially two sets of equations must 

be solved. The first relates the density- 
matrix components of the atomic system to 

the properties of the incident radiation 
field and the second describes the properties 
of the emitted radiation in terms of the 
density matrix. An inversion procedure, 
incorporating a transformation from the coordi¬ 
nate system used conventionally in the theory 
to that of the Instrument, is then used to 
determine corresponding magnetic field values. 
The Bommier code simplifies the computational 
procedure by calculating the field for a 
height of 50 arc sec and scaling this value 

for any given observational height, whereas 
the Landi Deg 1 1 Innocenti code computes £ 
in increments of 10 arc sec in height, and 
interpolates between heights. The Landi 
Degl' Innocenti code also incorporates an 
interpretation of Stokes V-profile data, 
which consists of an independent intensity¬ 
like component, d>(X ), and a Zeeman component, 
6$/6X (Landi Degl' Innocenti, 1982). This 
interpretation proves to be cruical in sifting 
out a single value for the magnetic field 
from a typical set of solutions, as discussed 
in Section 4 below. Since the Bommier code 
' was ; available in a fully-developed form, 
it'was used almost exclusively for the results 
reported here. 


4. Results 

4.1. Multiple Solution and Sifting Procedures 

It has been verified that synthetic 
linear polarization data, used as input, 
result In unique magnetic field vector solu¬ 
tion, or at most a symmetrical pair of solu¬ 
tions. The ambiguity is intrinsic to scatte¬ 
ring polarization in the presence of a magnetic 
field (see Landi Degl 1 Innocenti, 1982). 
However, when observational data are used 
multiple-valued solutions are generally produced 
that satisfy the fundamental requirement 
that a vector magnetic field predicted by 
one component of the D 3 profile is also predic¬ 
ted by the other. The inversion procedure 
of the code searches for the best solutions 
based on this criterion, in terms of a so- 
called distance function, D, where 

D = /^[ (A Q .) + (A U i ) 2 ] 

and where Q i = (Q^heo - Q^bs)/!, the summa¬ 
tion index running on the two components of 
D 3 line. The lack of uniqueness apparently 

result from residual noise, spatial and 
temporal averaging and imprecision In the 
geometrical parameters, all inherent in the 
observational data. 

Various procedures have been used to 
sift out Invalid solutions. The most straight 
forward is where the V-profile Is well defined. 
Theoretical V-profi les can be constructed 
from the linear polarization solutions and 
compared with the observed profiles. This 
procedure pan be sufficient to isolate the 
correct solution. As pointed out in the 
previous section, the V-profile consists 
of a component proportional to intensity 
and a Zeeman component, proportional to the 
derivative of the intensity, the two components 
denoted here by V and V 1 , respectively. 
The presence of V 1 therefore is evidence 
of a longitudinal component . of- the field 
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and Us absence that the field Is very nearly 
transverse. Further, the signs of Q, U, 
V and V' can be compared for self-consistency 
for a given set of B, h and x> and many solu¬ 
tions rejected on this basis. As an extension 
of this procedure, polarization maps for 
Q, U, V and V', for both components, so-called 
altitude diagrams, have been prepared. These 
are contour maps of polarization on a 6, 4 
plane for a given B, h, x set * intercom¬ 
paring a set of corresponding maps, inadmiss¬ 
ible solutions can be readily identified. 

By such methods, it is frequently found 
that all tentative solutions can be eliminated 
except a symmetrical pair. The correct direc¬ 
tion cap then be identified, at least with 
modest confidence, by comparison with photo- 
spheric magnetic fields. Uncertainty arises 
primarily from the possible substantial evolu¬ 
tion of a field distribution between the 
time of the prominence observations and that 
of the magnetograph record. Finally, self- 
consistency provides a further, general check 
on the validity of a set of solutions represen¬ 
ting measurements in the same prominence, 
especially if such measurements are repeated 
on two or more days. 

4.2. Summary of Magnetic Field Solutions 

The prominence measurements discussed 
here constitute a particularly useful set 
in that the profiles are generally of good 
quality (36 of 44 of the observed points 
are interpretable in the linear polarization) 
and extend ovbr four consecutive days, during 
which time the prominence gradually evolved 
in form and systematically increased in height. 
Since this prominence was relatively high 
and extensive and was observed over a four- 
day period, estimates of h and x could be 
obtained with more confidence than for typical 
prominences. 


The complete set of magnetic field 
vector solutions are presented in Figure 
5. The dashed outline of the prominence, 
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Fig.5. Vector field orientations and positions 
in the 19-22 April, 1970 quiescent proml- 
nonce. The ruimber sets are Meld strength, 
B (gauss), the .polar angle, 0 (degrees), 
measured from the local radius, and 
the azimuthal angle, <j>(degrees), around 
the local radius. Zero azimuth points 
west, +90° points north. 

imprecise because of the changing form of 
the prominence over the four-day period, 
Is intended to give simply a general frame 
of reference. Note the tendency for large 
values of |B| to occur in the leg of the 
prominence and that most of the remaining 
above-average values occur in the central 
part of the bulk of the prominence. The 
mean value is 23.8 G. 

Figure 6 is a latitude-longitude plot 
of field vectors for the same data set of 
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Fig.6. L atltude-longltude plot of field vectors 
for the 19-22 April, 1978 prominence, 
Arrows are proportional to the field 
strength. Surface magnetograph polarities 
are shown. 
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Figure 5. The length of the arrows are propor¬ 
tional to the field strength. The dashed 
line represents an approximate outline of 
the prominence as seen as a filament. The 
plus and minus signs indicate the surface 
magnetograph polarities. With some marked 
exceptions, the tendency for the horizontal 
component of the field vector to He close 
to 90° from the long axis ,of the prominence 
is obvious. 

5. Conments 

It has been pointed out above that 
the polarimeter field aperture, 3.5 arc sec 
by 10 arc sec, is large compared to the fine 
threads of emitting material, ~1 arc sec or 
smaller, that are sometimes observed in high 
resolution images of prominences. If such 

threads represent corresponding fine-scale 
magnetic structure, rather than simply 
being evidence of thermal instabilities, 
the observed mean polarization would be 
larger than that representative of the magnetic 
elements and the magnitude of the deduced 
magnetic field values would be correspondingly 
reduced. Implicit in the procedure for magnetic 
field derivation then is the assumption that 
the magnetic field is homogeneous over 
the polarimeter angular aperture. But even 

if the magnetic field is uniform over the 
field sampled by the polarimeter, line of 
sight contributions are likely to have a 
mixture of field components, which would 
result in a diminution of observed polarization 
and an increase to an unknown degree of inferred 
field strengths. The further assumption 
that the prominence material is effectively 
stationary during the period of observation 
is also not necessarily valid, and if so, 
the effect is also to increase the inferred 
field strengths. 


Conclusion 

Stokes polarimeter spectral profiles 
of He I Dg, representing many observational 
points within a quiescent prominence, have 
been Interpreted using a detailed quantum- 
mechanical theory of resonance polarization 
for the multlplet. The multiple set of 
solutions obtained when the linear polarization 
information alone is used can be reduced, 
generally to a symmetric pair, when the full, 
generalized theory is applied that interprets 
also the weak circular polarization Stokes 
profile. In most cases a unique solution 
can be found through additional information 
available from disk field polarity. 

Most of the derived magnetic field 
vectors are close to horizontal and near 
normal to the prominence surface, with some 
extreme exceptions, thus generally supporting 
the Kippenhahn-Schluter (1957) model. The 
values obtained for the magnetic field vectors 
are thought to be mostly accurate representa¬ 
tions of the actual prominence fields, a 
mean for this prominence of 24 G. 
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DISCUSSION: 

BHATTACHARYYA: (I) How do you calculate 

the strength of the magnetic field? Are 
there any calibrations from laboratory measure¬ 
ments? 

SMARTT: (1) The strength of the magnetic 

field is calculated by applying the theory 


of the Hanle effect for He I 0 3 line to the 
observed polarization. There has been no 
laboratory verification. 

BHATTACHARYYA; (2) Cannot one measure the 

field in prominences by Zeeman method? 

SMARTT: (2) Earlier attempts to measure 

magnetic field strengths in prominences based 
on the Zeeman effect, one now recognized 
as being invalid due to the dominance of 
the Hanle effect in the magnetic field regime 
appropriate for prominences, and especially 
so for lines which are not optically thin 
in prominences. 

VINOD KRISHAN: Have there been any attempts 
to produce analytic functional forms of the 
magnetic field from the derived numbers? 

SMARTT: So far, we have obtained only the 

vector magnetic field functions. We plan 
to investigate statistical relationships 
of the cumulative data and especially to 
use this information to try to understand 
the observed morphology of quiescent prominences. 
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MODELLING OF PLASMA STRUCTURES IN THE SOLAR ATMOSPHERE 

V, Krishan 
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Abiitrunt 

A maximal hallclty principle Is proposed to study the stability of plasma structures 
In tho solar atmosphere. This principle derives from the statistical treatment ol the magno- 
tohydrotJynomical turbulence. Tho steady state of tho plasma Bystem results from^ subjecting 
the system to the Invariance of total energy, tho magnetic heliclty and toroidal and pololdal 
magnetic fluxes. The statistical distributions of the velocity and magnetic fields are 
formulated. 


1. Introduction 

High resolution optical, UV, CUV and 
X-ray observations have established the 
spatial complexity of the solar corona. 
Active regions appear in the form of loops 
or arches. The spatial structure of these 
loops outlines the magnetic field geometry. 
From the observations of lines of Ca II, 
HI, CII, ' CIII, CIV, Ne VII, Mg VIII, Mg 
IX and Mg X, it is concluded that the loops 
consist of a cool core surrounded by a hot 
sheath which merges with the hot corona, 
Vaiana and Rosner (1978). 

2. Maximal Helicity Principle for 
Stability 

We assume plasma to be in the form 
of a cylinder of radius R and length L with 
rigid perfectly conducting walls at r = R. 
The MHO equations for an incompressible medium 
are; 

P[|f + (v.V) V ] * Is (VxB) X i-Vt-P9 (•!) 


ff * V x (VxB), V.V * V.B « 0 (2) 

where, the symbols have their usual meaning. 
The boundary conditions on B and V are; 

V r (r«R) « 0 « B r (r a R) (3) 

The geometry is indicated In Fig.l. 

Using the Incompressibility condition in 
the divergence of equation (1), one gets 
a relation between p, V and B as; 



Fig.l. Tho georriotry of a aolnt corona loop, 
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V 2 p = 1/4 tt V.[(VxB)XB] 


V-[p(V.V)V] + V.(pg) (4) 


which determines the density and temperature 
profile of the plasma. The three invariants 


of the plasma system are: 

Total energy E - /d 3 X (V 2 + B 2 ) (5) 

Magnetic he.11 city H M = /d 3 X A.B (6) 

Cross helicity H c = /d 3 X V.B ( 7 ) 

The poloidal flux ^ A 

r 7 0 5 2 

= constant ab r = R (8) 

•f 

where A is the vector potential, 


The velocity and magnetic fields are expanded 
in terrns of Chadrasekhar-Kendall functions 
given as: 


IK' 


a(n,m,q) = • '[ 4 * + ^ I «n.M> 

6 mk 

* e e [ -if ' - * (n ' m,c|) (9) 

* *,i i «■«> 


where 

■ W, r) e« + V 

A { n,m, q ) = ± CY* m£i + 

K ! n»0i tli 2,... 

n L 

m=0, ±1, ±2,... 

Y_ corresponds to that solution which 

makes a satisfy the boundary condition 
a r = 0 at r = R. This gives for m 2 + n 2 > 0. 


Now, all the fields and the invariants can 
be expressed in terms of a by using 

B * £ Cnmq A(n,m,q) £(n,m,q) a(n,m,q) 

( 12 ) 

V = £ Cnmq \(n,m,q) n(n,m,q) a(n,m,q) 

n,m,q 

(13) 

The statistical mechanics of the V and B 
fields is then formulated in a phase space 
whose coordinates are the real and imaginary 
parts of £'S and n‘S. The most probable 
state of the system is then constructed subject 
to the constraints of conservation laws. 
This leads to a canonical ensemble which 
has an exponential term having a lagrange 
multiplier for every conserved quantity. 
Lagrange multipliers are determined by requi¬ 
ring a given value of the expectation value 
of the conserved quantity, Montgomery et 
al. (1978). One chooses the canonical distri¬ 
bution subject to the invariance of E, H M , 
H c , vp t and ^ as s 

D * Constant exp [-aE - y’H c - 6^] (14) 

where a, 3, y and 6 are the Lagrange multi¬ 
pliers. For initially quiescent conditions 
H c =0 and, therefore y' =0. One can factor 
out the probability distribution for each 
expansion coefficient. The probability distri¬ 
bution for q(n,m,q) is: 

Pr)(n,m,q) = Kn,m,q exp [ ~a\ 2 (n,m,q) 

|n( n » m »q)| 2 ] ( 15 ) 


Rk y J‘ (y R) ■ + mX(n,m,q) Jm(Y R)=0 

n nmq m v 1 nmq ' ' x 1 nmq 

• ; ••• (3o) 

The eigen values for m = n = 0 are determined 
from the ratio of toroidal to poloidal flux: 


*t . -R 
% ' T 


Y 01 (y R ) 

qoq O ooq 

no,o,q) Jo'(y ofj(j R)" 



From this we find 

< jn(n.m.q)| 2 > * a" 1 X" 2 (n,m,q) (16) 

Similarly, one can carry out the procedure 
for other coordinates ?(n,m,q). One notices 
from equation (16) that even for initially 
quiescent conditions there is a finite mean. 
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square velocity for every mode in the system. 
Four scalar quantities, <E> , < H^> , < t > 
and < t|» p > are sufficient to determine 
the level of excitation of every mode In 
the system, 

3. Application to Coronal Loops 
3.1. Radial Temperature Variation 

We propose that the steady state of 
the loop can be described in terms of a single 
Chandrasekhar-Kenda 11 function character¬ 
ized by m = n = 0. Equation (4) for this choice 
can be solved to find the radial variation 
of the plasma temperature as; 

m. 

T{r) * -i X- (0,0,l)|n 2 (0,0,1)| C 2 (0,0,1) X 

k B 

C x { J 0 Z (x) + Jf (x) } - J,(x) Jo(x)] + To and 


3.2. Statistical Mechanics of V and B Fields 

In order to formulate the statistical 

mechanics, we include two values of the 

parameter q which labels the various roots 

of equation (11) as X Q , X v etc. One can 

express the average values of E, H M and 

in terms of X'S and the Lagrange multipliers. 

The Lagrange multipliers can be solved for 

by assigning some given values to <E> 

<H M >, and <ip > . As a numerical 

m t p g 

example, say L = 5x10 cm, iK/tfi = 1 

t p 

then equation (11) gives X Q « 2.2/R 

and X ^ - -2.6/R. The values of a , 6 and 

6 are determined from the following relation¬ 
ships; 

« = [zu-Xjiy (E-Xjiyr 1 x 


ST 

3F 


X s (0,0,1) In 2 (0,0, 1)1 c 2 (0,0,i) x 


J*(x) 

X 


(17) 


where T = T Q at r »0* x =Xr, XIs determined from 
equation (11). One can substitute for T(r) 
from eq. (17) in the line flux calculations 
and determine the spatial widths of the lines 
ClI, CIII, OIV, OVI, Ne VIII and Mg X which 
have mean temperatures of formation 5xl0 4 K, 
9 x10 4 K, 2x10 5 K, 3x10 5 K, 6x 10 5 K and 1.4xl0 6 K res¬ 
pectively. Now X 2 n 2 (0,0,1) is a measure of 
the turbulent velocity in the fluid andean 
be estimated from the observed Doppler widths 
of these lines. For X 2 n 2 ~10^(cm/sec)^ and 
T Q »-6 x 10 3 K, X = 2/R, R .» 10 9 cm, B ~1.5 Gauss, 
we find the spatial widths to be 0.02 X"V 
0.025 X" 1 , 0.12 X’ 1 , 0,1215X*’ 1 , 0.3X _1 and 
> 0.5X" 1 for CII, CIII, OIV, OVI, Ne VII and. 
Mg X respectively, This behaviour certainly 
conforms to the cool core and hot sheath type 
of loops. •' 


C - ( - 3E + | (X 0 + Xj) H m } ± 

{ E 2 -|(X 0 + X 1 ) 2 H- 2 - BX^r 2 


- E (X 0 + Xj) H m l 1 / 2 ] 

X, 


(18) 


8 = 


and 


1 


2(E - 1/a - Xjiy 2{E - 1/a - X 0 H ) 


, . . r «i (E - l/« - X x H m ) 

' ** [ w*i>— 

ti (E - 1/a - X 0 H„,) , 

X^ X 1 - Xq) 

where 6 0 = 2 ttRC(0,0,1) Y 001 Ji(y 0Q1 R) 

5l = 2 RC(0, 0> 2) Yqd2 Y 002 R) 


E, 

for 


(19) 


( 20 ) 


H m and.. jjj'j, are the observables, Since 
the case of coronal loops, there are 
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no direct measurements of these quantities, 
we shall have to fix their valies from other consider 
rations. - For example i can be obtained 
from energy balance arguments which give 
E -40*8 ergs (Levine and Wlthbroe 1977). 
We choose ■ ~ Q,Q5xl0 3? ergs cm, The two 

values of a are - 2 . 04 x 10 "^ (ergs)" 1 ar >d 
0.98x10"^ (ergs)” 1 . We find a ~ 2.04x10^ 
(ergs)to be appropriate for the approxima¬ 
tions involved in arriving at (18). Eq. 
(19) gives p ? 0.2xl0" 3 ^ (ergs cm)' 1 , and 

6 = 10 ^ (Gauss cm 2 )” 1 . A representative 

choice of ij> t ~10 17 Maxwell gives 6 = 10" 19 
Thus the distribution function for this 
particular Indicative example is : 

D = Const e” 2 ' 04E/E x 

e -0.01^ t / iji t (21) 


We observe that total energy has a rather 
narrow distribution compared to the distribu¬ 
tions of the magnetic heliclty and toroidal 
flux. More details about the application 
to coronal loops are presented in the paper 
by Krishan (1983a,b). 
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Abstract 

The solar dectimotric radio emission 'consists of fl variety of fine structures. This 
reflects the complexity of the corona ot this level. Two observations, viz. absorption 
bursts and short duration narrow bond bursts nro presented. An interpretation of each 
of those observations Is provided. 


1. Introduction 

The solar decametrlc radiation emanates 
from the region of 1-2 solar radii. The 
decametrlc noise storms consist of a continuum 
which may last for one to several hours. 
Superimposed over this continuum are a variety 
of fine structures. The most often observed 
features are the type III, type I lib, drift 
pairs, S bursts and so on (Sawant 1982). 
However, there are large periods of time 
when the continuum alone is observed (Gergely 
and Kundu 1975). In addition to the above 
mentioned fine structures, Sastry et a 1. (1983) 
observed sudden reductions in the decametric 
continuum called the absorption bursts. 
Also,, short duration narrow band bursts which 
have a small drift rate were observed recently 
(Gopalswamy et al. 1983b), In this paper, 
we discuss the absorption bursts and the 
short duration narrow band bursts and provide 
a physical interpretation for them. These 
observations provide Information about the 


complex magnetic structure in the outer corona. 
The decametric corona is one region to which 
sufficient attention has not been paid. 

The significance of this region Is reflected 
In the following facts* (i) the co-rotation 
of the corona ceases at this level (11) the 
magnetic field configuration changes predomi¬ 
nantly from closed to open field structures 
(iii) new spectral features appear which 
are absent in the higher frequency region. 

Our observations do not always have 
flare associations. Hence, the physics of 

these observations pertains .to the'quiescent' 
corona by which we mean that the changes 
in magnetic field are small and may correspond 
to the signatures of an evolving active, region, 
The examples we are considering seem to ref lect 
the ‘transition’ character' of the magnetic 
field in the decametric corona. The short 
duration narrow band bursts which occur re Tati - 
vely rarely • may corresponds,to the existence 
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of closed magnetic field configurations while 
the absorption bursts Indicate open field 
lines. Both these configurations are possible 
as it is a transition region. 

2. Absorption Bursts 
2.1. Observations 

The absorption bursts are seen as dips 
in the continuum flux for a short duration. 
Detailed observations are provided elsewhere 
(Sastry et al. 1983, Gopalswamy et al. 1983a), 
The most important characteristics can be 
summarised as fol lows: 

(i) The absorption duration, on an average 
is about 2-3s, 

(ii) The absorption bursts occur over a relative 
bandwidth of about 5%, 

(iii) The absorption bursts do not have any 
drift in the frequency time plane, 

(iv) Large number of absorption bursts have 
a simple time structure. However, there 
are double and triple structures observed 
and sometimes multiple time structures, 


(v) The intensity reductions are > 40$, 

Typical examples are shown In Fig.l. The 
frequency observed is around 34.5 MHz. 

2.2. Interpretation 

The sudden reductions are also observed 
in high frequency levels, (e.g. Benz and 
Kuijpers 1976, Fokker 1982). The cause of 
absorption depends upon both the generation 
and propagation characteristics of the conti¬ 
nuum radiation. Since the radiation generation 
and their propagation outwards are different 
at different frequency regions, the interpreta¬ 
tion should take account of these facts. 
So, one can classify the interpretations 
into two catagories: 

(i) Absorption at the generation stage 
where the process of plasma wave 
generation or the process of conversion 
of plasma wave into electromagnetic 
radiation could be stopped temporari ly 
to cause the observed reduction, 
(ii) Absorption could be attributed to 
the propagation effect in that a 
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Fig.l. Eight channel record of absorption burst obtained on April 13, 1981. The 
ordinate is the intensity of radio emission. The absorption bursts are at 
069323 UT and 065209 UT. 
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'screen' Is present on the path of 
. the radiation for a short .while that 

causes the intensity reduction. 

Benz and Kuijpers (1976) regarded 
that the sudden reduction in the decimetric 
continuum must be due to filling of loss 
cone of electron distribution trapped in 
the closed magnetic configuration which caused 
the decimetric continuum. While the loss 

cone is filled by the newly injected electrons, 
the loss cone instability for plasma waves 
is stopped and hence one sees a reduction 
in the radiation. Fokker (1982) pointed 
out the drawback of Benz and Kuijper's theory 
and proposed that the reduction could be 
due to inhomogeneities created on the path 
of radiation by solitons or shocks which 
impinge on the flux tubes Into which the 
continuum radiation is ducted. The problem 
with Fokker's Idea Is that it needs open 
flux lines at decimetric levels which are 
quite rare. 

If one wants to apply these theories 
to decametric continuum, one faces many diffi¬ 
culties. The decametric continuum is explained 
as a result of the Rayleigh scattering of 
enhanced plasma fluctuations caused by a 
diffuse electron beam along open magnetic 
field lines (Levin 1982). Since the theory 
of Benz and Kuijpers needs a closed field 
configuration, the loss-cone filling mechanism 
is not appropriate here. Fokker's (1982) 
idea contains the correct ingredients but 
one has to be sure of the ducting of continuum 
radiation at decametric region. There is 
no concrete evidence for: the existence of 
low density plasma In the coronal streamers 
to cause the radiation "ducting. Moreover, 
the inhomogeneities created by shocks■ ire 
solitons of small scale lengths which may 
not be able to reduce the flux by S0%. Adher¬ 
ing to the idea of Fokker (1982) one must 
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consider an Important aspect of the shock 
waves such as the creation of low frequency 
turbulence in the shock front an<^ its wake 
(Galeev, 1976). The extent over which the 
plasma turbulence exists is sufficiently 
large (~ 10^ cm) so that the large decametric 
continuum source of size 10' (Kundu, 1983) 
could be screened -temporarily. When the 
continuum radiation 'passes through this region 
of turbulence, it interacts with low frequency 
turbulence (ion-sound turbulence is the one 
generated in the shock front and exists in 
the wake with a high levelllf energy density). 
The result is the generation of Langmuir 
waves at the expense of the continuum radiation. 
Once the Langmuir waves build up to a certain 
level, the reverse interaction of generation 
of radiation, becomes., important and hence 
the saturation occurs in 'the absorption. 
The ion-sound turbulence decays bec'ause there 
Is no free energy available for its generation 
after the shock passage. Once the intensity 
of ion-sound turbulence falls below a critical 
level, there is no effective conversion of 
continuum radiation, into,.Langmuir waves and 
hence the intensity recovers to the original 
continuum level. This critical ..level of 
ion-sound ; turbulence could be estimated from 
the requirement .that the optical depth of 
the turbulence layer exceeds unity. From 
weak turbulence theory one can obtain the 
optical depth of the turbulence layer and 
calculate the critical level In terms of 
the effective temperature T of the ion- 

S q -i, 

sound turbulence; T $ > 2.4xKrK. From the 
observed bandwidth and the scale length of 
coronal inhomogeneity, one can obtain the 
linear extent of the turbulence layer as 
~10 cm. The collisiona 1 damping time of 
the ion-sound turbulence is in the right 
range of observed duration of absorptions. 

'The non-drif ting-character of-the absorp- 
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tions support the fact that shock waves are 
propagating perpendicular to the radial magne¬ 
tic field and hence perpendicular to 
the line of sight, A schematic model is 
proposed in F1 g-(2>, for the origin of such 



The short-duration narrow band bursts 
are identified as type I bursts at low frequen¬ 
cies. The solar origin of these bursts is 
established In several ways. Interference 
due to lightning, static ignition, etc. are 
of broad band with zero frequency drift. 
'.The broadcast interference is of narrow band 
but again, it has no frequency drift. The 
characteristics of- these bursts could be 
summarised as follows: (see Fig.3).' 


24 July 1981 34-50MHZ 



" (•/ 

Fig. 2. A modal lor Iho absorption process. : 

a shock wave. When an emerging magnetic 
field line opens in a neighbouring site, 
it can drive a low Mach number shock across 
the open field lines. The multiple -fine 
structures could be interpreted as the passage 
of successive shocks. The different time 
structures would then correspond to the differ¬ 
ent time intervals between successive shocks. 
The mathematical details have been given 
by Gopalswamy et a 1. (1983a). 

If the shock propagates through the 
source region itself, then the enhanced plasma 
fluctuations will get converted into radiation 
which propagates perpendicular to the line 
of sight and hence one can observe a reduction 
in the continuum flux. 

3. Short-duration Marrow-band Bursts 
3.1. Observations 


34-30 MHz 



Fig. 7. Eight channel record obtained on 

July 24 1901. The ordinate Is the 

Intensity qI radio emission. 

(i) The frequency drift is very small 
which 1$-250 KHZ S" 1 . 

(ii) The single frequency duration is 
I -300 ms. 

(ii1) The bandwidth is. very narrow ~3 %. 
(iv) The single frequency duration increases 
with decreasing frequency. The total 
increase is -100 ms, 

(v) The bursts are superimposed over weak 
continuum. 
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The short duration narrow band bursts 
cannot be S-bursts (McConnel 1982) as the 
latter have very large drift rates. Though 
the drift rates are slntilar to type 1 ' II' bursts, 
one can rule out the possibility of type 
II because the type II 1 s have a very large 
single frequency duration of several seconds. 
These bursts therefore have to be type I 
bursts. Since the occurence of closed magnetic 
fields in decametric region is relatively 
rare, the frequency of occurence of the type 
I bursts is also rare, which is consistent 
with the observations. 

3.2. Interpretation 

The bursts could not have been generated 
by electron cyclotron mechanism which is 
acceptable only when the plasma frequency 
is less than or close to the gyrofrequency 
of the electrons. The burst generation due 
to electron beam could also be ruled out 
because in this case the drift rates will 
be extremely large and contradict the observa¬ 
tions. There is another possibility of getting 
the drift rate*, a radio signal generated 
at a particular level is of finite extent 
in time, So, different Fourier components 
of the signal travel with different velocities 
because the refractive index is different 
for different frequencies. Fo‘r such a process 
to be effective, a calculation of the dispersion 
measure shows that, the electron column density 
is extremely small and one may need six orders 
of magnitude more. Based on these facts 
one can clearly come to the conclusion that 
the bursts must be created by a weak shock 
wave moving out in the corona and can explain 
the observed slow drift. 

Spicer et al. (1981) proposed the emerging 
flux theory to explain the , type I bursts 
and the same thing is regarded as a working 
model (Wentzel 1982). The basic ingredients 


of the theory are two wave modes and a hoi 
particle distribution. The high frequency 
wave mode is created by the hot particles. 
The low frequency mode produces the hot parti* 
cle distribution by stochastic acceleration. 
The interaction between the two wave, modes 
produces the radiation, ' The scenario is 
as follows. Newly emerging magnetic field 
pushes against the pre-existing coronal itiagne- 
tic fields and produces shocks. In the shock 
front and its wake, a high level of lower 
hybrid (Lll) waves are generated. Part of 
these LH waves stochastically accelerate 
electrons which are subsequently trapped 
in the closed magnetic arches and develop 
loss-cone anisotropy. The free energy In 
this electron distribution is fed into the 
upper hybrid waves (UH). ,A weak shock of 
slightly super Alfvenic velocity can generate 
the condition necessary for radio-emission 
and the drift rate of the bursts is the velo¬ 
city of the shocks. The emerging flux theory 
can explain our observations. If the observed 
bursts are caused by a weak shock then the 
density jump across the shock should be direct¬ 
ly related to the observed bandwidth. Therefore 
the observed bandwidth and the estimated 
shock velocity (from the drift rate) could 
be used in conjunction with the Rankine-Hugon- 
lot relations to estimate the ambient magnetic 
field at decametric levels. The result' is 
of the order of 1 G which is in very good 
agreement with other estimates such as frequen¬ 
cy splitting and polarization measurements 
of type I bursts (ter Haar and Tsytovlch 
1980). This agreement is an indirect confirma¬ 
tion of the emerging flux theory in addition 
to some observational evidence provided by 
Karlicky’and Giricka (1982). 

According to the emerging flux model, 
type I source is very tiny and the nablation 
is emitted instantaneously (Wentzel 1982). 
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The observed duration must be due to the 
finite extent of the generated radiation 
beam In which the rays making larger angle 
with the direction of the magnetic field 
lag behind, those at smaller angles* This 
mechanism explains the short duration of 
the bursts. As the shock proceeds to lower 
frequencies, the extent of radiation beam 
may increase and also, there could be an 
increase in the range of angles in which 
the radiation is emitted. These two changes 
can cause increase in duration with decrease 
in frequency as observed. The weak continuum 
over which the bursts are superimposed must 
have been created by the interaction of UH 
waves with some unknown low frequency fluctua¬ 
tions. 

In fact a weak shock can generate 

a variety of low and high frequency turbulences 
depending upon the strength of the shock 

and the coronal conditions which the shock 
faces as it propagates through the corona. 
The high frequency turbulence (Buneman Instabi¬ 
lity) is excited only at the Initial portion 
of. the shock front where the gradients are 
steep, Among low frequency waves,’ the LH 
waves are generated for isothermal conditions. 
The ion-sound waves are generated only when 
there is non-isothermality. Because of this 
the generation of ion-sound turbulence is 

always viewed with suspicion. But if there 
is an initial steep gradient., then the Buneman 
Instability generated can be quenched by 
essential heating of electrons and thereby 
causing non-isothermality (Tidman and Krai 1, 
1971). In fact, the characteristic time 

scale for electron heating through Buneman 
instabi lity is t b - /3(m e /m|)^ 3 w" 3 - 10' 7 s. 
This is very small compared to the shock 

transit time (« pe .V^/C)"* ~5.8 X 1 6 5 s (Vj - 
shock velocity). Hence the non-isothermality 
will be met within about 1/100th of the shock 


thickness. Under such a situation, therefore, 
ion-sound Instability, could be excited, 
A detailed calculation of the ion-sound insta¬ 
bility has been reported elsewhere (Gopalswamy 
et al. 1983b). 

For the UH waves to provide the adequate 
brightness temperature, the low frequency 
waves should have an energy density (Melrose, 
1980); 

W q 6/3 V e c 

nT e ~ 71 %e L n V ph 

where • ; 

W q energy density of low frequency 
nT e mode o interacting with the UH 

waves, 

n - electron density of the corona 

■ V e - electron thermal velocity 

V ph - phase velocity of the low frequency 
waves 

t n - coronal density gradient scale 

length 

c - velocity of light . 

For n = 10 8 cm" 3 and T g = 10 6 K, the above 

inequality becomes ; 


nT 


> 2.4 x 10 


-6 


The energy density of LH waves Is calculated 
by Spicer et al. (1981) as ; 


W LH 


(Jj 


ce 

‘pe 


AB 

B 


'A 


-1 




where V. is Alfven velocity, AB is the relative 

' - b 

jump in magnetic field across the shock, 






Finn Structure In Baler Dflanmatrlo Radio Emission 


53 


Is the Alfv^nic Mach number, (m Q /m^) is 
the ratio of electron to Ion mass and C. 

5 

Is the sound velocity. Following the similar 
procedure one can calculate the energy density 
of Ion-sound waves by extending the one-dimen¬ 
sional calculation of Krall and Book (1969) 
to three dimensions (Gopalswatny et al. 1983b) 
asi 

Wl „ C S V A AB JVi 2 A, 1/2 

nT e B Mft Z-B 'me' 

x [1 + In - 1 . 

UT 

ce 

For identical conditions of the shock In 
the decametrlc corona, 

^ “ 7.7 x 10' 6 and ^ « 2.8 x 10" <1 , 

e e 

so that 


This estimate shows that the LH waves only 
marginally satisfy the energy density require¬ 
ment while the ion-sound wave density is 
much larger than the threshold value. Moreover, 
the Ion-sound waves are generated over a 
wide range of wave numbers compared to that 
of LI! waves, thereby satisfying the resonance 
condition more easily, 1«e M the overlap 
between UH waves will be more in the case 
of Ion-sound waves. Then, subject to the 
condition that the ion-sound waves are genera¬ 
ted , they are better candidates compared to 
the I.H waves. 

Another important aspect of the emerging 
flux theory Is the generation of UH waves. 


In the theory of Spicer ■ et. al. (1981), the 
IH waves stochastically accelerate the elect¬ 
rons parallel to the magnetic field. . But 
UH waves are generated by a loss-cone distribu¬ 
tion of electrons. Therefore, the accelerated 
electrons, travelling essentially parallel 
to the magnetic field should convert a major 
part of their parallel energy into perpendicu¬ 
lar energy. Vlahos et al. (1982) have proposed 
anomalous Doppler resonance instability 
to generate a loss-cone distribution through 
pitch angle scattering. However, if the 
low frequency turbulence involved is due 
to ion-sound waves, there is another possibi¬ 
lity of generating high frequency plasma 
turbulence. Since the ion-sound frequencies 
are of the same range as the whistlers, they 
can get non-linearly scattered into whistler 
waves in a time-scale of 0.01s (Kaplan & 
Tsytovich, 1973). The whistler polarization 
is such that the energy of the electrons 
is increased predominantly in the direction 
perpendicular to the magnetic field due to 
the whistler absorption. The time-scale 
over which the transverse energy of the electi¬ 
on increases is of the order of 0.5 s (Kaplan 
& Tsytovich, 1973). Since both these time 
scales are much less than the collisions! 
clamping time of the Ion-sound waves, an aniso¬ 
tropic electron distribution could be produced, 
the free energy content of which could be 
fed into high frequency plasma oscillations 
through a mechanism similar to one dlsgussed 
by Levin (1982). 

In conclusion we would like to point 
out that while the emerging flux theory can 
adequately explain our observations as well 
as other type I phenomena, the low frequency 
models that provide the opacity could be 
different depending upon the exact conditions 
prevailing in the shocks and the corona. 
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4. Conclusions 

We hatfe presented some observations 
of the fine structures in the decametric 
solar radio emission. We have provided repre¬ 
sentative events which reflect the complex 
structure of the decametric corona. Our 
observations and their interpretations indicate 
the importance of this special region of 
the corona in understanding the radio physics 
of the Sun. As the decametric corona has 

not received much attention from theorists, 
a detal led investigation of the radio processes 
from this level will provide an opportunity 
for a better insight into outer corona as 
well as interplanetary phenomena. 
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DISCUSSION: 

BHATTACHARYYA: In all the cases mentioned 

for absorption events did you find any suitable 
loop configuration from optical pictures? 

GOPALSWAMY: Yes, that is what I was just 

discussing with Dr. R.N.Smartt, But, we 
have not specifically identified any event. 
He has promised to send some photographs 
in this regard, I am planning to search 
from our own optical records, 
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CONCLUDING DISCUSSION 


R. Howard 
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The subject of this Colloquium is 
appropriate for the 75th anniversary of the 
discovery of the Evershed effect. It is 
also a timely and important subject. Today 
we have heard contributions on a variety 
of topics, and we have seen that there remain 
a number of important unsolved questions. 
We have learnt a great deal about the Inter¬ 
action of plasma and magnetic fields In recent 
decades, but the phenomena of the solar atmo¬ 
sphere continue to defy physical explanation 
and their field remain as challenging as 
the day 75 years ago when Evershed discovered 
motions in sunspots. 

Mathematical modelling of the small 
scale magnetic structures of the solar atmo¬ 
sphere is leading to more accurate physical 
descriptions of such features, but at the 
same time the refinements in such analyses, 
as we have seen here today, bring with them 
the necessity for more detailed, high resolution 
observations of the solar surface. The promise 
of a quantum jump in spatial resolution from 
space observations lies still in the distant 
future, so far now the best efforts of ground- 
based observers will be needed to aid the 
theoreticians. Today we have heard from 
both theoreticians and observers, and we 

" .Pre«enlly ill National Solar Obnotvutory, Tucson, 
Arizona, U.S.A. Operwlud by the Association 
of Unlvflrnilltts lor Rnsenrch In Agronomy, Inc. 
under conltnct with the Notional Science Foundation. 


see that progress es reported from both groups, 
Meeting such as this which bring together 
these researchers are valuable means of advan¬ 
cing the field. 

The nature of the solar activity cycle 
remains one of the most puzzling questions 
in solar physics. Today, we have heard a 
fine review of this field and an interesting 
new model has been proposed. Few areas of 
astronomy have seen such a long and heated 
debate as that of the solar dynamo. In recent 
years the stellar astronomers have joined 
the 'fray with observations of stellar cycles. 
This has added a valuable new dimension to 
the subject. Today we saw observations relating 
to the cycle and one is tempted to suggest 
that although new observations are always 
desirable, in this field we may have too 
many observations. That is, the interaction 
has the difficult job of relating those large 
scale, long term phenomena that are truly 
relevant to the fundamental cycle mechanism 
and discarding those that are only incidental 
to the important physics of the process, 
It may be that we have already In hand the 
key to the solution of the cycle, if we only 
knew what it was and don't mean to suggest 
that we should stop observing large scale 
phenomena. Clearly the search must continue 
unti T we can make further progress. Currently 
the kinematic dynamo models are favoured 
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by many who have worked on them. But, as 
we have seen today, they have their weaknesses 
and many people feel that we must look elsewhere 
for a cycle mechanism. The search continues. 

It is appropriate that at this meeting 
we have heard a good review about the Evershed 
effect. We seem to have made little progress 
in understanding this phenomena in 75 years, 
and after all this time we are still faced 
with observations that seem contradictory. 
Of course perhaps not all sunspots show the 
same dynamical properties, but if this is 
the case, then we can add a further puzzling 
dimension to the Evershed effect. The Very 
highest possible spatial resolution is required 
to study the motions in sunspots. This is 
another field where a breakthrough may be 
anticipated when sub arc-second observations 
from space become available. In this case, 
as In others, we need high resolution velocity, 
magnetic field, and brightness information 
with extended coverage in time (at least 
some hours). Until such space observations 


become possible, and in fact, in preparation 
for them, we should use modern ground-based 
techniques to define better the dynamical 
properties of sunspots. In particular, it 
seems to me, differences in the amplitude 
or form of the Evershed effect should be 
Investigated. Are such irregularities related 
to the size, shape or age of the spots? What 
form does the Evershed effect take for irre¬ 
gularly shaped penumbrae? The answers to 
questions such as these can be obtained from 
the ground and should precede space observa¬ 
tions. 

We work in an exciting field that 
does not lack for Important unsolved questions. 
Among these are the nature of small scale 
motions in sunspots and elsewhere on the 
solar surface and the fundamental mechanism 
of the solar cycle. An important start in 
the study of these processes was made 75 
years ago here by Evershed, and the opportuni¬ 
ties for significant advances in our knowledge 
of these topics is even greater today than 
it was then. 





